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Detector  active  area 
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Adaptive  threshold 
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Eq.  (C9) 
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Background  (quantum  noise)  limited  photodetector 
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Background-structure-induced  fluctuations  in  the 
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Background-structure  limited  (performance  regime) 

C 

Eq.  (74b) 

Ratio  of  noise  variance  to  scene  radiance 

CFAR 
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d 
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Eq.  (C7) 
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False-alarm  rate 

v 


FPA 

Focal-plane  array 

FT 

Fixed  threshold 

hit) 

Eq.  (11) 

Impulse  response  of  the  postdetector  electrical  filter 

h(l) 

Eq.  (17) 

Time  derivative  of  h(t) 
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Eq.  (11) 

Transfer  function  of  the  postdetector  filter; 

Fourier  transform  of  h(i) 

ho  ( t ) 

Eq.  (31) 

Impulse  response  of  the  threshold  filter 
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Eq.  (B20) 

Fourier  transform  of  /»0(f) 

1, 

Eq.  (4) 

Subinterval  of  the  time  interval  /0 

lo 

Eq.  (3) 

An  interval  of  time 

1RST 

Infrared  search  and  track 

J 

Eq.  (2) 

Number  of  threshold  crossings  (an  integer  random  variable) 

K 

Thresncla-gain  constant 

*0 

Eq.  (A4) 

Ratio  of  the  focal-plane  irradiance  to  the  scene  radiance 

Eq.  (Cl) 

Expected  value  of  A  {i) 

mu 

Eq.  (Cl) 

Expected  value  of  B(t ) 

Eq.  (A4) 

Expected  value  of  the  focal-plane  irradiance 

mj 

Eq.  (2) 

Expected  value  of  the  crossing  count 

rhj 

Eq.  (3) 

Threshold-crossing  rate  (the  time  derivative  of  m, ) 

""o 

Eq.  (36) 

Rate  of  positive-slope  zero  crossings 

>»L 

Eq.  (74a) 

Scene  spatial  radiance  distribution 

mj  (f) 

Eq.  (A8) 

Two-dimensional  spatial  Fourier  transform  of  mL( r) 
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Eq  (A2) 

Mean  photon  flux  incident  on  a  photodetector 

mx 

Eq.  (9) 

Expected  value  of  the  unfiltered  current  X(i) 

»h 

Eq.  (12) 

Expected  value  of  the  video  signal  )'(/) 

VI 


m) 

Eq.  (13) 

Expected  value  of  Y(t ) 

m  f 

Eq.  (147b) 

Time-advanced  and  time-delayed  versions  of  mY(t) 

my 

Eq.  (60) 

Estimate  for  mY 

m,hY 

Eq.  (61) 

Expected  value  of  thy 

m  ■ 

"  > 

Eq.  (B29) 

Expected  value  of  &  y 

M 

Eq.  (151) 

Noise  bandwidth  parameter  of  the  postdetector 
electrical  filter 

MTF 

Eq.  (A6) 

Modulation  transfer  function 

N 

Eq.  (58) 

Number  of  series  TD1  detectors 

NEI 

Noise-equivalent  irradiance 

Pit) 

Eq.  (A7) 

"Pupil  function"  of  a  focal-plane  detector 

9(0 

Eq.  (A7) 

Two-dimensional  spatial  Fourier  transform  of  Pit) 
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Eq.  (6) 

Probability  of  detection 

P,a 

Eq.  (7) 

False-alarm  probability 

QNL 

Quantum-noise  limited  (performance  regime) 
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Two-dimensional  focal-plane  location  vector 
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Eq.  (A5) 

Photodetector  current  responsivity  [A  W'1] 
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PASSIVE  INFRARED  SURVEILLANCE, 
PART  II:  THRESHOLD-CROSSING  RECEIVERS 


INTRODUCTION 

The  objective  of  an  infrared  surveillance  receiver  is  to  passively  detect  targets  at  long 
range  while  keeping  false  alarms  acceptably  infrequent.  In  this  report  we  begin  with  a  review  of 
the  surveillance-receiver  analysis  model  developed  in  Part  1  [1],  which  relates  expected  system 
performance  to  appropriate  descriptions  of  the  background  scene,  the  electro-optic  sensor,  and 
the  structure  of  the  signal  processor.  We  then  apply  the  principles  elucidated  in  Part  I  to  the 
processor  design  problem.  Some  of  the  ideas  underlying  the  approach  and  some  of  the  main 
results  are  summarized  below. 

The  signal  processor  must  examine  the  preamplifier  output  waveform  and  classify  it  either 
as  targetlike  or  backgroundlike  in  structure,  according  to  a  set  of  criteria  chosen  by  the  system 
designer.  Since  the  structure  of  the  time  variations  in  the  video  signal  figures  so  importantly  in 
how  well  the  processor  performs,  we  now  discuss  the  mechanisms  that  give  rise  to  these  time 
variations. 

The  video  waveform  at  the  output  of  the  postdeteclor  filter  of  a  scanning  background- 
limited  (BLIP)  photodetector  displays  fluctuations  having  two  fundamentally  different  origins: 

•  Quantum  noise.  Time-of-arrival  fluctuations  of  the  individual  photons  incident  on 
the  detector  cause  random  variations  in  the  photosignal  called  quantum  noise  or  photon  noise. 

•  Background-structure-induced  fluctuations.  Scanning  the  sensor  across  spatial  varia¬ 
tions  in  the  scene  radiance  induces  time  variations  in  the  photosignal  referred  to  as 
background-structure-induced  (BS!)  fluctuations. 

Each  of  these  two  basic  mechanisms  enters  into  the  statistical  representation  of  the  video  signal 
in  complementary,  but  distinctly  different,  ways. 

The  background-structure-induced  fluctuations  enter  into  the  statistical  description  of  the 
video  signal  as  a  modulation  imposed  on  the  mean  value  of  the  video.  Photon  noise  enters  into 
the  statistical  representation  of  the  video  as  the  origin  of  the  variance  of  the  video  about  its 
mean.  The  statistical  characterization  is  complete  once  the  time-varying  mean  and  variance  of 
the  video  signal  are  specified,  since  the  video  is  a  nonstationary  Gaussian  random  process.* 


Manuscript  submitted  September  14,  1979. 

‘This  fundamental  description  of  the  video-signal  statistics  has  apparently  not  been  previously  applied  to  the  passive- 
surveillance  problem.  Validity  of  the  nonslalionary  Gaussian  model  for  the  video  is  assured  whenever  the  expected 
number  of  photons  incident  on  the  detector  during  a  sensor  dwell  time  is  much  greater  than  unity  (21. 
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Each  of  the  candidate  processor  designs  we  describe  in  this  report  consists  of  two  parts: 
one  part  for  estimating  tne  mean  video  signal  (as  determined  by  the  background  structure),  and 
a  second  part  for  estimating  the  noise  variance  due  to  quantum  noise.  A  main  result  in  this 
report  is  the  proposal  of  a  new  signal-processor  structure  for  generating  high-confidence  esti¬ 
mates  of  the  time-varying  rms  quantum  noise.  Implementation  of  the  proposed  quantum-noise 
estimator  requires  a  sensor  with  either  true  dc  response  or  synthetic  dc  restoration. 

The  second  main  component  of  the  signal  processor,  the  mean-video  estimator,  has  a 
more  critical  impact  on  system  performance  than  the  quantum-noise  estimator  when  operating 
against  highly  structured  scenes.  Unfortunately,  the  mean-video  estimator  is  also  more  difficult 
to  design,  since  the  designer  is  required  to  make  assumptions  regarding  the  nature  of  the  struc¬ 
tured  background. 


We  concentrate  our  discucrion  of  mean-video  estimators  on  three  types  of  linear  filters, 
evaluating  their  performance  against  a  succession  of  increasingly  complex  model  backgrounds: 

•  uniform  scenes  of  known  radiance, 

•  uniform  scenes  of  unknown  radiance, 

•  slowly  varying  nonuniform  scenes,  and 

•  scenes  containing  objects  of  arbitrary  contrast  and  size.* 

Improved  mean-video  estimators  can  no  doubt  be  developed  by  incorporating  actual  meas¬ 
ured  background  scenes  into  the  processor-design-and-evaluation  procedure  and  by  expanding 
the  scope  of  candidate  processor  designs  to  include  nonlinear  algorithms  as  well  as  linear  ones 
(31.  Although  a  representative  database  of  distributions  of  measured  background  infrared  radi¬ 
ance  is  thus  required  both  for  designing  and  assessing  more  advanced  processor  algorithms,  we 
do  not  consider  the  availability  of  such  data  in  this  report. 

In  addition  to  the  limited  scope  of  our  mode!  backgrounds  and  mean-video-estimation 
algorithms,  a  number  of  implicit  assumptions  underlie  the  approach  to  receiver  design  outlined 
above.  Perhaps  foremost  of  these  is  the  basically  deterministic  nature  of  our  model  back¬ 
grounds.  For  example,  the  only  requirement  on  each  member  of  the  class  of  "slowly  varying" 
structured  backgrounds  is  that  the  relative  change  in  radiance  from  one  sensor  dwell  time  to  the 
next  be  much  smaller  than  unity.  We  make  no  attempt  to  give  a  statistical  description,  such  as 
a  probability  density  or  correlation  function,  either  to  individual  scenes  or  to  the  ensemble  of 
scenes. 

The  utility  of  the  "unknown  but  nonrandom"  infrared-background  representation  is  ration¬ 
alized  by  paraphrasing  Van  Trees  [4,  p.  456]  (with  the  words  in  brackets  modifying  the  original 
text): 

"It  is  sometimes  unrealistic  to  consider  the  signal  ...  as  a  random  waveform. 

I'or  example,  we  may  know  that  each  lime  a  [false  alarm)  occurs  the  [back¬ 
ground  scene!  will  have  certain  distinctive  features.  If  the  (background)  is 

‘This  final  category  of  model  backgrounds  subsumes  the  class  of  real  backgrounds  containing  point  objects  (such  as  tar¬ 
gets)  of  cither  positive  or  negative  contrast. 
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modeled  as  a  sample  function  of  a  tandoin  process  we  may  average  out  the 
features  that  are  important  Here  it  is  more  useful  to  model  the  {back¬ 
ground]  as  an  unknown,  but  nonrandom,  waveform." 


The  spatial  correlation  function  (or  its  Fourier  transform,  the  Wiener  spectrum)  is  an 
example  of  a  statistical  scene  description  that  illustrates  this  point.  Typically,  the  presence  or 
absence  of  a  point  target  in  a  given  scene  has  a  dominant  effect  on  system  perf  •  mance  while 
having  a  negligible  effect  on  the  correlation  function  of  a  large,  structured  scene 

Since  we  require  a  deterministic  description  for  the  background  spatial  variations,  we 
sometimes  refer  to  oi*  performance  model  and  our  performance  assessments  as  background 
conditional.  Another  implicit  assumption  in  our  design  procedure  concerns  the  restricted 
nature  of  the  target/background  discriminant  that  we  use. 

Of  all  the  possmle  target/background  discriminants,  such  as  radiance  variations  with 
time*,  space,  polarization,  and  spectral  content,  we  will  consider  only  the  spatial  discriminant. 
Moreover  we  will  make  no  attempt  to  exploit  the  spatial  character  of  backgrounds.  Any 
attempt  to  do  so  is  complicated  by  the  enormous  diversity  of  structure  displayed  by  natural 
scenes.  Instead  we  will  key  on  the  known  point  character  of  the  target.  Each  time  we  propose  a 
candidate  processing  structure,  we  will  ask  the  same  question:.  Are  there  any  nonpointlike  struc¬ 
tural  features  that  could  induce  a  threshold  crossing?  If  so,  the  signal  processor  will  be 
modified  in  a  way  intended  to  correct  the  apparent  deficiency.  Our  ultimate  objective  is  to  con¬ 
struct  a  receiver  that  is  sensitive  only  to  the  presence  of  pointlike  objects  in  the  sensor’s  field  of 
view. 


At  present  the  limiting  efTect  on  system  performance  is  generally  background  structure 
rather  than  d:tector  sensitivity.  Thus  improvements  in  detector  sensitivity  are  of  little  use, 
unless  we  can  trade  off  sensitivity  for  improved  performance  against  structured  scenes. 

In  the  final  part  of  this  report  we  will  examine  the  usual  search-set  rule  of  thumb  relating 
electrical  bandwidth  and  sensor  dwell  time.  We  will  find  that  the  use  of  targe  video 
bandwidths,  in  violation  of  the  rule  of  thumb,  enables  us  to  sacrifice  some  degree  of  perfor- 
man: .  against  uniform  backgrounds  in  the  interest  of  enhanced  performance  against  nonuni¬ 
form  backgrounds.  Fortunately  the  performance  degradation  thus  incurred  against  uniform 
scenes  can  be  offset  by  means  of  focal-plane  detector  arrays  (FPA)  employing  time-delay-and- 
integration  (TDD  logic.  Thus,  we  will  establish  a  concrete  link  between  advancements  in  com¬ 
ponent  technology  and  overall  system  performance, 

BACKGROUND 

Concepts  and  Definitions 

The  two  conventional  parameters  for  measuring  how  well  the  infrared  search-and-track 
(1RST)  sensor  accomplishes  its  dual  objectives  of  target  detection  and  false-target  rejection  are 
the  probability  of  target  detection  (P0)  and  the  false-aiarm  rate  (FAR) 


'For  example  tracking  algornnms  are  neither  modeled  no;  dixcusseil,  even  thougn  tliev  ma>  pla>  an  important  rule  in 
largel/background  discrimination  as  well  as  target  tracking  i5| 
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The  1RST  performance  model  described  below  [1]  requires  as  input  complete  descriptions 
of  the  IRST  sensor  and  the  scene  radiance  distribution  (possibly  including  a  target).  As  output 
the  model  generates  the  probability  that  the  IRST  device  declares  (rightly  or  wrongly)  a  target’s 
presence  in  the  scene.  A  priori  knowledge  as  to  whether  a  target  was  in  fact  present  in  the 
specified  scene  allows  interpretation  of  the  target  declaration  probability  (my)  either  as  a  target- 
detection  probability  ( PD )  or  as  a  false-alarm  probability  (PFA). 

Photon  fluctuation  noise  is  the  only  stochastic  aspect  of  the  model;  the  background  must 
be  specified  as  a  radiance  map  of  arbitrary,  but  deterministic  structure.  Thus  the  performance 
predictions  made  with  this  model  are  background  conditional. 

The  IRST  analysis  model  that  follows  applies  to  threshold-comparison  receivers  configured 
as  in  Fig.  1..  The  current  X(t)  at  the  output  of  the  detector  is  input  to  an  electrical  filter  of 
transfer  function  //(/)..  The  output  current  Y(t)  of  the  electrical  filter  is  compared  with  a 
threshold  Tof*)--  If  TO)  exceeds  the  threshold,  the  presence  of  a  target  is  declared;  otherwise 
no  target  declaration  is  made. 


Fig  1  -  A  basic  ihreshoid-comparison  receiver  The  photodeiector 
in  this  figure  is  idealized  in  that  it  is  presumed  to  have  a  perfect  all¬ 
pass  electrical  fr. quern. y  characteristic,  the  frequency-dependent  part 
of  the  detector  responsitivity  is  lumped  together  with  the  transfer 
function  of  the  postdetector  filter  to  obtain  H(J)  A  target  declara¬ 
tion  is  made  whenever  the  filtered  current  Yit)  exceeds  the  thres¬ 
hold  level  y0(r) 


Figure  1  suggests  that  the  threshold-establishing  mechanism  should  suppress  clutter- 
induced  threshold  crossings  (false  alarms)  by  increasing  T0O)  when  TO)  is  clutterlike.  Rather 
than  allow  T0O)  to  take  on  an  a  prion  constant  or  functional  value,  it  is  necessary  to  establish 
the  threshold  by  some  means  that  adapts  T0O)  to  the  prevailing  background  conditions. 

Proper  operation  of  the  processor  requires  that  Mr)  <  T0O)  when  there  is  no  target  in 
the  sensor’s  field  of  view  A  threshold-crossing  event  is  defined  as  taking  place  at  any  time  in 
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such  that  y(rrc)  —  y0(rrc).  AH  threshold  crossing  detectors  are  presumed  to  be  sensitive  only 
to  positive-slope  crossings,  for  which* 


Hire)  **  Yo 


(1) 


Since  the  current  Y(t)  is  a  random  process  due  to  time-of-arrival  fluctuations  of  the  light 
quanta  incident  on  the  detector, 

•  threshold  crossings  are  random  events, 

•  the  threshold-crossing  times  tTC  are  random  variables,  and 

•  the  number  of  threshold  crossings  ( J)  that  occur  during  a  specified  time  interval  :s 
an  integer  random  variable. 


We  define  the  average  value  of  J  as  the  mean  crossing  count  my 

mj  -  E[J)( 


(2) 


where  the  statistical  expectation  operator  F{)  denotes  an  average  value  over  the  photon 
fluctuation  statistics  of  the  incident  light  [1]. 

We  now  introduce  a  threshold-crossing  rate  m,  to  help  explain  the  relationship  between 
mj  and  the  conventional  search-sct  performance  measures  PD  and  FAR  [1],  By  definition 


mj  &  J  mj(t )  di. 


(3) 


where  /0  is  an  appropriately  chosen  time  interval. 


It  is  generally  found  that  the  principal  contributions  to  the  crossing-count  integral  in  Eq 
(3)  accrue  in  the  neighborhood  of  a  discrete  set  of  times  which  in  turn  can  each  be  associated 
with  a  structural  feature  in  the  scene,  As  the  IRST  sensor  scans  across  a  target  (or  a  back¬ 
ground  feature  that  simulates  a  target),  the  crossing  rate  mj(t)  becomes  a  sharply  peaked, 
pulselike  waveform. 

We  now  define  the  time  interval  I,  as  tne  entire  period  of  time  during  which  the  sensor  is 
scanned  across  a  single  structural  feature  in  the  scene  (such  as  a  local  maximum  in  the  radiaice 
distribution).  By  presumption,  /,  is  a  subinterval  of  the  total  scan  interval  /0  in  Eq.  (3).  It  can 
be  shown  that  the  increment  in  associated  witn  the  interval  /,,  that  is, 

Am,(/,)  =  f  rhj(t)  di.  (4) 


is  generally  a  number  be-ween  0  and  1: 


0  <  Am//,)  <  1. 


(5) 


We  show  in  Part  I  111  that  A m//,)  may  be  interpreted  as  the  probability  of  a  thieshold  crossing 
sometime  during  the  interval  /,,  that  is,  as  the  probability  of  a  target  declaration  during  the  time 
interval 


*B>  convention  a  dotted  variable  denotes  the  time  derivative  ol  me  corresponding  undolled  variable 
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If  the  structural  feature  scanned  during  the  interval  /,  is  known  a  priori  to  be  a  target,  we 
interpret  the  probability  of  target  declaration  as  a  probability  of  target  detection: 

PD  -  Am;(/,).;  (6) 

Conversely,  if  the  IRST  sensor  is  known  a  priori  to  be  scanning  a  region  in  space  that  does  not 
contain  a  target  during  we  interpret  the  quantity  Am//,)  as  the  probability  of  a  false  alarm 
curing  the  interval  /,:. 

PFAU,)  -  Am//,),  (7) 

Dividing  the  total  scan  interval  /0  into  a  succession  of  disjoint  subintervals  /,,  each  of  which 
brackets  an  identifiable  structural  feature  in  the  scene,  we  can  write  the  overall  false-a'arm  rate 
as  the  expected  number  of  false  alarms  during  /0  divided  by  the  duration  of  /0: 

FAR  -  j±—  £  />,//,).  (8a) 

Thus,  from  Eq.  (7), 

FAR  -  -  -  £A m//,)  (false  alarms  per  second) ,  (8b) 

*  * *0'  | 

where  H/o)  is  the  duration  of  ‘ime  interval  /0  (seconds)  and  PfA  (/,)  -  A m//,)  is  the  probabil¬ 
ity  of  a  false  alarm  during  time  interval  /,, 

In  summary,  target  detections  and  false  alarms  are  both  manifested  by  the  signal  proces¬ 
sor  as  threshold  crossings.  Consequently  we  can  express  the  search-set  performance  measures. 
PD  and  FAR  in  terms  of  a  threshold-crossing-rate  function  m/r). 

Following  the  next  subsection  we  will  formulate  the  crossing  rate  m;  in  terms  of  appropri¬ 
ate  characterizations  of  the  IRST  sensor  and  the  background  scene  under  observation.  We  will 
develop  the  formulation  for  the  threshold-crossing  rate  mj  in  two  steps  111: 

•  First,  we  will  express  the  quantities  m/r)  and  as  functions  of  the  scene  radiance 
distribution  and  various  optical  and  electrical  attributes  of  the  IRST  device.  We  will  define  the 
quantities  m/r)  and  c ry(i)  as  the  ensemble  mean  value  and  variance  of  the  filtered  current 
K(f)  (Fig.  1),  Analagous  expressions  for  the  mean  and  variance  of  Y{i),  Fo(r),  and  Y0(t)  are 
also  required,  and  we  will  present  them.  All  such  mean  values  and  variances,  individually  and 
collectively,  are  referred  to  as  the  current  statistics  of  the  IRST  sensor. 

•  We  will  express  m/r )  in  terms  of  the  current  statistics  (as  a  function  of  m)(o>,  etc.). 

The  expression  for  the  threshold-crossing  statistic  m,  in  terms  of  the  statistics  of  Y(t)  and 
K0(r)  is  devoid  of  physical  content,  it  applies  equally  well  to  any  nonstationary  Gaussian  pro¬ 
cess.  We  will  introduce  all  of  the  physical  parameters  -  the  optical  and  electrical  characteristics 
of  the  IRST  device  and  the  radiance  of  the  background  scene  -  through  the  formulation  for  the 
current  statistics. 

Current  Statistics 

The  primary  current  statistic  of  interest  is  the  average  value  m/r)  of  the  detected  photo¬ 
current  A(r)  (Fig.  1)' 

m/r)  =  E[X(t)\.  (9) 
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As  in  Appendix  A,  mxit)  can  be  written  as  a  function  of  the  scene  radiance  distribution  and  a 
number  of  important  IRST  sensor  parameters. 

The  detected  current  Xit)  and  filtered  current  Yit)  are  related  by  the  linear  system- 
input/output  relation 

Yit)-  f"  h(t-T)X(r)dT  =  h(t)®XU).-  (10) 

J  —oo 

The  symbol  ®  is  a  shorthand  notation  for  the  convolution  process  defined  by  the  integral 
expression  in  Eq.  (10),  Also,  hit)  in  Eq.  (10)  is  the  impulse  response  of  the  postdetector  filter; 
that  is,  hit)  is  the  Fourier  inverse  of  the  filter  characteristic  Hif)  in  Fig.  1: 

Hif)  -  hit)  exp  i- jin  ft)  dt.  (11) 


When  the  expected  value  is  taken  of  both  sides  of  Eq.  (10),  it  follows  directly  from  Eq. 
(9)  that 


myit)  =  £(y(r))  -  hit)®mxit). 

(12) 

The  following  relationship  can  also  be  found  directly  from  Eq.  (10): 

myit)  =  E{Yit )}  -  myit).. 

(13) 

Expressions  similar  to  Eq.  (12)  are  also  required  for  the  quantities 
where  by  definition 

< ryit )  and  <r2it), 

(r2it)  =  £T{ly(/)— mK(f)]2) 

(14) 

and 

<Xy(r)  =  EUhO-m^r)]2}. 

(IS) 

An  adaptation  of  Eq.  (4.3.13)  in  Ref,  6  leads  to 

<Tyit)  -  eh2it)®mxit), 

(16) 

where  e  is  the  electronic  charge  and  hit)  and  mxit)  are  the  same  quantities  as  in  Eq.  (12).  As 
in  Appendix  B  of  Part  I  (1 J,  it  follows  from  Eqs.  (10)  and  (16)  that 

cr  >•  ( r )  -  e[h  (r)l2  ®  mxit),  (17) 


For  the  special  case  of  spatially  uniform  scenes,  the  analysis  of  Appendix  A  shows  that 
mx  is  time  independent.  It  then  follows  from  Eq.  (12)  that 

mY  ”  mx f  hit)  dt.  (18) 

However,  from  Eq.  (11) 

HiO)  -  f  hit)  dt  (19) 

-oc 

Thus,  from  Eqs.  (18)  and  ('9) 


mt  “  Hi0)mx 


(20) 
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for  uniform  backgrounds.  For  bandpass  transfer  functions,  H( 0)  -  0.  It  follows  from  Eq.  (20) 
that  Y(t)  is  a  zero-mean  process  for  the  special  case  of  uniform  backgrounds  and  bandpass 
postdetector  titters. 

Similarly  it  can  be  shown  from  the  Parseval  relation  for  Fourier  transforms 

J*~  hHt)  dt  -  J_“  \H(/)\1  df  (21) 

that  Eqs.  (16)  and  (17)  simplify  as  follows  in  the  uniform  background  limit: 

i emx  J  m  \H(f)\2  df  (22) 

and 

cr  y  *■  emxf  m  (2 irf)1  |// (/)  |2  df.  (23) 

Assuming  that  the  transfer  characteristic  //(/)  has  been  normalized  to  a  peak  value  of 
unity,  that  is, 

max  //(/)-  1,  (24) 

we  define  the  conventional  noise  bandwidth  of  //(/)  as 

A/  =  Jo“  \H(f)\2  df.  (25) 

It  follows  from  Eqs.  (22)  and  (25)  that 

«r  y  “  2 emx  A/,  (26) 

which  is  the  usual  Poisson  shot-noise  variance  formula  (given  in,  for  example,  Ref,  7),  By 

inference  Eq.  (16)  is  a  generalization  of  the  usual  shot-noise  variance  expression  to  nonstation¬ 
ary  shot-noise  processes. 

The  motivation  for  the  results  presented  in  this  subsection  is  provided  in  the  following 
subsection,  in  which  we  express  the  threshold-crossing  rate  riij(t)  as  a  function  of  the  current 
statistics.  (It  should  be  recalled  from  the  Concepts  and  Definitions  subsection  that  the  usual 
search-set  performance  measures  PD  and  FAR  may  then  be  calculated  in  terms  of  the  crossing 
rate  rhj.) 

In  summary,  the  current  statistics  mr,  mY,  crr,  and  crr  may  all  be  written  as  relatively 
simple  functions  of  the  postdetector  filter  impulse  response  hit)  and  the  average  value  mx(t ) 
of  the  unfiltered  current  X(t)  (Fig.  1): 


mr(t)  =  h  (r)  ®mx(t), 

(27) 

mY(t)  =  mY{t), 

(28) 

a-yit)  =  eh2(t)  ®mx(t), 

(29) 

o-2y(t)  =  e[h(t)]2  ®mx(t). 

(30) 

Expressions  for  the  threshold-current  statistics  can  be  written  directly  by  analogy  to  Eqs.  (27) 
through  (30)  for  the  threshold  approach  depicted  in  Fig.  2: 
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mYoU)  -  Kh0U  -  Td)  ®mxU), 

(31) 

mY()U)  -  mYoU). 

(32) 

cr  yQ(f)  -  eK  h${t  -  Td)  ®mx(t), 

(33) 

<r2y  U)  -  eK  [h0  (t  -  Td)]2  ®mx(t), 

(34) 

where  K  and  Td  are  the  gain  and  time  delay  respectively  in  the  threshold  circuitry  and  /i0O)  is 
the  Fourier  inverse  of  the  threshold  filter  characteristic  H0(f)  (Fig.  2), 


Fig.  2  —  A  simple  adaptive-threshold  scheme.  The  transfer  function 
txp(~j2nf  Td)  introduces  a  delay  of  Td  seconds,  ensuring  decorrelation  of  the  ran¬ 
dom  processes  Y(t)  and  Y0(t)  Figure  l  shows  the  significance  of  the  random  sig¬ 
nals  XU),  YU),  and  Y0U). 


The  mean  unfiltered  current  mx(t)  that  figures  so  prominently  in  Eqs.  (27)  through  (34) 
is  expressed  in  Appendix  A  as  a  function  of  the  scene  radiance  distribution  and  a  number  of 
important  IRST  sensor  parameters. 

Threshold-Crossing  Rates 

We  will  now  briefly  outline  the  method  developed  in  Part  I  [11  for  calculating  the  average 
number  of  times  mj  that  the  random  current  YU)  crosses  the  random  threshold  Y0(t)  during  a 
specified  time  interval  (Figs.  1  and  2).  The  relationship  between  the  mean  crossing  count  mj 
and  the  usual  search-set  parameters  P0  and  FAR  was  briefly  discussed  in  the  Concepts  and 
Definitions  subsection  of  this  report  and  discussed  in  greater  depth  in  Part  1  (1). 

We  first  express  the  crossing  count  as  the  time  integral  of  a  threshold  crossing  rate  rhjU) 
(Eq.  3).  We  then  write  the  crossing  rate  m,  as  a  function  o?  the  current  statistics  as* 

m/t)  -  mJoU)  exp[— w2(r)/2l .  (35) 

where  by  definition 

mJoU)  =  (2rr)_l  (o-  YYJcr  yYi}),  (36) 


•Equation  (35)  for  my  is  not  strictly  valid  in  the  neighborhood  o''  "mean-crossing  times"  (to  be  defined  by  Eq  (42))  A 
more  general  formulation  for  mj  valid  everywhere  is  derived  in  Part  1  111  However,  as  we  will  note  following  Eq 
(42),  the  more  general  expression  for  is  not  usually  required,  since  it  integrates  to  unity  in  the  neighborhood  of  a 
mean  crossing  [l,2l 
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b 

(37) 

2  —  2,2 
^no  ~  arY  +  <TY0> 

(38) 

u(t)  =  ( my0  -  my)/(Tyy0. 

(39) 

As  indicated  by  Eqs.  (4),  (6),  and  (8),  evaluation  of  the  important  search-set  parameters 
PD  and  FAR  requires  the  time  integration  of  Eq,  (35)  to  obtain  the  mean  crossing  count  mj 
over  proscribed  intervals  of  time: 

mj  -  J*  rhj(t)  dt,  (40) 

where  1\  is  the  time  interval  of  interest. 

When  mxU)  is  slowly  varying  compared  to  /;(/),  it  can  be  shown  that  »?,(/)  is  practically 
constant  with  time  (as  will  be  expressed  by,  Eq.  (48)).  In  this  case,  it  follows  from  Eqs  (35) 
and  (40)  that 

mj  *  mJ()  J  exp[-u}(/)/2]  dt,  (41) 


As  we  will  now  discuss,  the  integration  in  Eq.  (41)  is  often  quite  simple  to  perform. 


Mean-Crossing  Times  and  Times  of  Closest  Approach 

Experience  with  the  numerical  evaluation  of  Eq,  (41)  has  shown  that  the  principal  contri¬ 
butions  to  the  integral  accrue  in  the  neighborhood  of  a  discrete  set  of  times.  Moreover,  it  has 
been  shown  that  these  important  discrete  limes  are  of  two  types,  "mean-crossing  times”  and 
"closest-approach  times"  (1,2], 

As  illustrated  in  Fig.  3a,  the  mean-crossing  times  tmi  are  those  times  that  simultaneously 
satisfy  the  two  conditions 

"rOme)  -  rnYoUmc)  (42) 

and 

>  my0  Unc)..  (43) 

For  each  solution  of  Eq.  (42)  that  satisfies  constraint  (43),  that  is,  each  time  the  mean  current 
my(t)  crosses  the  mean  threshold  wl#(/)  with  positive  slope,  the  expected  number  of  crossings 
mj  is  incremented  by  unity.  Whenever  mean  crossings  exist  during  the  interval  /,,  it  is  gen¬ 
erally  not  necessary  to  perform  the  integration  in  Eq.  (40);  in  this  case  the  expected  number  of 
threshold  crossings  mj  is  well  approximated  by  the  number  of  mean-crossing  times  rm  during 
the  interval  f.  Thus,  the  expression  for  thj  is  generally  not  needed  to  evaluate  the  contribu¬ 
tions  to  mj  that  occur  in  the  neighborhood  of  mean  crossings,  since  each  mean  crossing  gives 
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rise  to  a  threshold  crossing  [1].*  Clearly  it  is  desirable  that  no  mean  crossings  occur  except 
when  there  is  a  point  target  in  the  scene;  this  may  be  taken  as  the  first  principle  of  search-set 
design  for  operation  against  structured  backgrounds. 


Fig.  3a  -  A  "mean-crossing  time"  denoted  tm .  A 
mean  crossing  occurs  when  the  mean  video  signal 
mY  crosses  the  mean  threshold  signa.  myQ  with 

positive  slope.  Asymptotic  evaluation  of  the 
crossing-rate  integral,  Eq.  (40),  shows  that  the  pro¬ 
bability  of  a  threshold  crossing  approaches  unity  in 
the  neighborhood  of  a  mean-crossing  time  Ill. 


If  mY(t)  lies  below  myQ(t)  on  the  time  interval  /1(  that  is,  if  there  are  no  mean  crossings 
during  /,,  the  crossing-count  integral,  Eq.  (40),  is  generally  dominated  by  contributions  accru¬ 
ing  in  the  neighborhood  of  "times  of  closest  approach.”  As  illustrated  in  Fig.  3b,  closest- 
approach  times  are  those  times  that  satisfy  the  equation 

mY(tca)  -  mY(tca).  (44) 


Fig,  3b  —  A  "clOEsst-approach  time,"  denoted  /l0 
The  threshold-crossing  rate  mj  typically  increases 
many  orders  of  magnitude  in  the  immediate  neigh-- 
borhood  of  a  closest -approach  time. 


Fixed-Threshold  Processors 

For  fixed-threshold  (FT)  processors 

cr  v-  -  cr  y  -  0 
'o  • o 


and 


">r0-y  0- 


(45a) 

(45b) 


"Although  a  threshold  crossing  is  assured  in  the  near  neighborhood  of  a  mean  crossing,  the  threshold-crossing  time 
tj('  is  not  equal  to  the  meati-crossing  time  lmi since  //•(■  is  a  random  variable  (as  staled  following  Eq  (1)1,  while  t,m 
is  deterministic 
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It  follows  from  Eqs.  (35)  through  (39)  and  (45)  that 

rhj(t)  -  (2ir)_1  (<ry/<Ty)  cxp[-w2(f)/2],  (46) 

where 

u(t)  -  (y0- mr)/(Ty.  (47) 

The  current  statistics  my,  cry,  and  cry  in  Eqs.  (46)  and  (47)  are  obtained  from  Eqs.  (27) 
through  (30),  and  y0  in  Eq.  (47)  is  a  fixed  threshold  value. 

When  mx(t)  is  slowly  varying  compared  to  h(t),  it  follows  from  Eqs.  (23)  and  (26)  that 

(«  °o  \  1/2 

(A/)-'Jo  f2\H(f)\2  dj\  =  (A/)rmj.  (48) 

Equations  (41)  and  (48),  taken  together,  provide  a  relatively  simple  formalism  for  calculating 
the  expected  crossing  count  of  a  fixed-threshold  receiver  operating  against  a  slowly  varying 
background  scene. 

In  the  limit  of  time-invariant  mx  (perfectly  uniform  backgrounds)  Eq.  (41)  simplifies  still 
further  to 

mj  -  (A f)rm  exp( — AT2/2)  <  r(/,).  (49) 

where  TU\)  is  the  duration  of  /,,  (A/),mj  is  once  again  given  by  Eq.  (48),  and  the  fixed  thres¬ 
hold  y0  has  presumably  been  adjusted  such  that 

u  -  (y0-  my) /ay  -  K,  (50) 

Equation  (49)  is  originally  due  to  Rice  [8)  and  has  often  been  applied  to  the  calculation  of 
false-alarm  rates  for  infrared  search-set  operation  against  uniform  scenes  [9,10], 


Initial  Assessment  of  Adaptive  Vs  Fixed-Threshold  Processors 

Some  of  the  advantages  and  disadvantages  of  adaptive-threshold  (AT)  processing  as  com¬ 
pared  with  fixed  threshold  (FT)  processing  may  be  evaluated  by  comparing  Eqs.  (39)  and  (47), 
appropriate  for  AT  and  FT  receivers  respectively. 

It  follows  from  Eq.  (47)  that  an  FT  processor  suffers  a  false  alarm  whenever  the  peak  tar¬ 
get  amplitude  is  less  than  the  peak  clutter  amplitude.  This  situation  is  depicted  in  Fig.  4a. 
Thus  target-to-clutter  ratios  less  than  unity  cause  the  FT  processor  performance  to  be 
"background-structure  limited"  (BSL).  In  this  case  each  false  alarm  can  be  associated  with  a 
structural  feature  in  the  background.  The  effect  of  quantum  noise  (as  reflected  in  the  magni¬ 
tude  of  tr  r,  for  example)  is  then  totally  overshadowed  by  background-structure  effects. 

Inspection  of  Eq.  (39)  and  Fig.  4b  shows  that  an  adaptive-threshold  (AT)  processor  need 
not  suffer  background-induced  mean  crossings.  If  the  filter  H0(f)  (Fig.  2)  can  be  chosen  such 
that  mYo(t)  tracks  the  background-induced  variations  in  w >  (/),  the  background-induced  mean 
crossings  can  be  eliminated.  Although  false  alarms  may  still  occur  as  random  events  due  to 
quantum  noise,  the  probability  is  generally  small  that  any  given  structural  feature  in  the  scene 
will  cause  a  false  alarm. 
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Fig.  4a  —  A  critical  shortcoming  of  fixed-threshold  processing.  The 
slowly  varying  maximum  centered  at  tp  has  its  origin  in  the  nonuni¬ 
form  background  scene  The  narrower  peak  at  /0  is  due  to  a  target. 
The  likelihood  of  a  false  alarm  grows  rapidly  as  the  threshold  level 
>'o  is  reduced,  There  is  no  way  for  tne  constant-threshold  processor 
to  detect  the  target  peak  at  r0  without  also  incurring  a  false  alarm 
arsing  from  the  clutter  peak  centered  at  ip.  A  plot  of  the 
threshold-crossing  rate  mjU)  corresponding  to  this  figure  would 
show  that  the  probability  of  a  threshold  crossing,  and  hence  a  false 
alarm,  is  far  greater  at  time  ip  than  at  any  other  time,  since  ip  is  a 
time  of  closest  approach. 


Fig.  4b  —  An  important  advantage  of  adaptive  threshold  pro¬ 
cessing  The  mean  current  m^t)  is  the  same  as  for  Fig,  4a. 
The  adaptive  threshold  }0(r)  accurately  tracks  the  slowly  vary-- 
ing  background  signal  but  not  the  more  rapidly  varying  target 
signal.  Thus,  target  detection  is  assured,  and  the  probability  of  a 
false  alarm  is  kept  acceptably  small.  As  contrasted  with  the 
situation  of  Fig  4a,  the  probaoility  of  a  false  alarm  is  now  no 
greater  in  the  neighborhood  of  time  ip  than  at  any  other  time. 


An  IRST  sensor  operating  in  this  regime  (such  as  in  Fig.  4b)  is  said  to  be  'quantum-noise 
limited"  (QNL)  in  its  performance..  In  this  context  QNL  operation  is  clearly  preferable  to  BSL 
(background-structure-limited)  operation. 
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Uniform  Scenes  of  Known  Brightness 

The  adaptive-threshold-performance  advantage  just  described  is  realized  only  when  the 
background  scene  is  nonuniform.  The  performance  of  an  adaptive -threshold  processor  is  infe¬ 
rior  to  the  performance  of  a  properly  designed  FT  processor  when  :he  background  scene  is  uni¬ 
form  and  of  known  brightness.  In  this  case  the  adaptive-threshold  fabe-alarm  rate  (FARat)  is 
greater  (worse)  than  the  fixed-threshold  false-alarm  rate  (FARpr): 

FARat  >  FARft..  (51) 

For  the  false-alarr.i-rate  comparison  expressed  by  Eq.  (51)  to  be  meaningful,  the  adaptive- 
threshold  gain  K  in  Fig.  2  has  been  adjusted  to  achieve  equal  target-detection  sensitivities  for 
the  two  processors  being  compared, 

We  now  define  uniform-background  false-alarm-penalty  (FAP): 

FAP  =  10  logic  (FARAT/FARpr)  dB,  (52) 


Uniform -Background  False- A larm -Penalty 

The  uniform-background  false-alarm  penalty  is  a  measure  of  the  performance  disadvan¬ 
tage  that  accrues  when  an  AT  processor  is  used  against  a  uniform-background  scene  of  known 
brightness.  As  shown  in  Part  1  (11,  the  false-alarm  penalty  of  the  processor  structure  in  Fig.  2 
may  be  written  as 

FAP  -  2.17 
where  by  definition 

a  =  A/q/A/  (54) 

and  where  A/ and  A/0  are  defined  as  the  noise-equivalent  bandwidths  of  Hif)  and  //0(/), 
respectively  lEq.  (25)  and  Fig.  2). 


mYo  -  mY 


(Ty 


+  51og|0(l  -  a(l  -  a)]  dB. 


(53) 


As  an  ,xample,  we  assume  that  the  threshold  constant  K  in  Fig.  2  is  adjusted  until  "the 
threshold  is  live  sigmas  above  the  mean,"  that  is. 


">y 

(Ty 


5. 


(55) 


when  the  search  set  is  observing  a  uniform  scene  of  known  brightness.  Furthermore  we 
assume  that  the  noise  bandwidth  of  the  target  filter  Hif)  in  Fig.  2  is  twice  as  large  as  the  noise 
bandwidth  of  the  threshold  filter  Hoif)'- 


a  “ 


A/q 

A  / 


1 

2  ‘ 


It  follows  from  Eqs.  (53)  through  (56)  that  the  false-alarm  penalty  is 


(56) 


FAP  -  17.46  dB, 


(57) 
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corresponding  to  a  value  of  FARat  (adaptive-threshold  false-alarm  rate)  about  56  times  worse 
than  FARpr  (fixed-threshold  false-alarm  rate).  Clearly  this  is  a  severe  degradation  in  perfor¬ 
mance. 


Time  Delay  and  Integration 


Fortunately  the  disadvantage  of  AT  processing  can  be  more  than  recovered  by  means  of 
time  delay  and  integration  (TDI).  The  effect  of  TDI  is  simply  to  multiply  the  mean  current 
mx(t)  by  the  number  N  of  series  TDI  detectors.  It  follows  from  Eqs.  (27)  and  (29)  that  TDI 
multiplies  my  and  <ry  by  the  factors  N  and  NUi  respectively.  It  can  be  shown  that 


FAPtd, 


1  +  a  —  N 
a 


2.17 


a 

1  +  a 


my0-  my 

<Ty 


+  51og,otl  -a(l  -a)]  (58) 


for  the  AT  processor  depicted  in  Fig.  2,  employing  N  series  TDI  detectors.  It  follows  from  Eqs. 
(55),  (56),  and  (58)  that 

FAPTD1  -  ((3  -  2A)  (18.08)  -  0.62]  dB  (59) 

for  the  example  given.  Thus  the  uniform-background  false-alarm  penalty  suffered  by  the  AT 
processor  of  Fig.  2  is  more  than  compensated  by  the  use  of  just  N  -  2  TDI  detectors. 
Apparently  the  false-alarm  rate  can  be  made  arbitrarily  small  by  increasing  the  number  of  TDI 
detectors.  Alternatively  the  false-alarm  rate  can  be  held  constant  with  increasing  N  by  making 
the  threshold  gain  K  (Fig.  2)  inversely  proportional  to  NU2..  The  TDI  performance  advantage  is 
then  realized  as  a  target-detection  sensitivity  enormously  increased  over  that  of  the  N  -  1 
single-detector  IRST  sensor.  IRST  performance  against  uniform  scenes  improves  monoKu- 
cally  as  the  number  of  TDI  elements  is  increased. 


Again,  the  false-alarm  penalty  is  meaningful  only  as  a  characterization  of  AT  processor 
performance  against  the  most  benign  type  of  background— a  uniform  background  of  known 
brightness.*  However,  the  motivation  for  employing  AT  processing  lies  in  the  unfortunate  * e?!- 
ity  that  natural-background  scenes  are  often  highly  structured. 

In  following  sections  of  this  report  we  will  characterize  a  number  of  cauidate  proofing 
schemes  with  respect  to  their  performance  against  a  succession  of  more  c.  nip'.  K  scenes.  We 
discuss  the  rationale  underlying  the  design  of  these  candidate  processors  m  the  following  sub¬ 
section. 


Ideal  (CFAR)  Adaptive-Threshold  Receivers 

The  candidate  processor  designs  described  in  this  report  each  consist  of  two  parts:  one 
pan  for  estimating  the  mean  video  signal  my,  and  a  second  part  for  estimating  the  rms  quan¬ 
tum  noise  cr Trier  signal  processor  establishes  the  adaptive  threshold  }'„(/)  as  shown  in  Fig.  5: 

r((i)-«,(nH(i,(i)(  (60) 


•The  relevance  ui  scene  raJiance  mt  to  search  set  performance  against  umlorm  scenes  follows  directly  from  bus  (26) 
and  (AS)  The  nuise  variance  </y  increases  linearly  with  scene  radiance  Thus  operation  against  bright  uniform 
scenes  occasions  decrees :d  target  detection  sensitivity  and/or  increased  talse  alarm  rale  as  compared  with  operation 
against  umlorm  scenes  of  lesser  brightness 
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Fig.  S  —  Generic  structure  of  the  signal  processors  considered  in  this 
report.  The  expected  value  and  rms  value  of  the  video  signal  Y(t) 
are  denoted  as  mYU)  and  try 0)  respectively.  The  signals  rhy(t)  and 
&YU)  are  estimate:  for  mY(r)  and  c ry(t).  The  adaptive-threshold 
constant  K  is  the  gain  of  an  ideal  amplifier.  As  shown  in  Fig.  1,  a 
target's  presence  is  declared  (rightly  or  wrongly)  whenever  the  video 
signal  Y(t)  exceeds  the  threshold  signal  Yq(i). 


where  mY  and  &Y  are  estimates  for  mY  and  cry  respectively..  The  adaptive,  threshold  constant  K 
in  Eq.  (60)  is  a  design  parameter.  The  estimates  mY  and  &Y  and  the  threshold  T0  are  all  sto- 


chastic  processes.  With  the  expected  values  of  mY  and  &Y  denoted  as 
tively,  it  follows  from  Eqs.  (39)  and  (60)  that 

m and  respec- 

u(t)  -  K(mit/<TYYo)  +  -  mY)/crYYo 

Equation  (61)  can  be  written  as 

(61) 

u(t)  -  K(  1  +  f-)  +ycv 

(62) 

\ 

where  the  normalized  mean  estimation  errors  and  «•  are  defined  as 

"1  y  U  y 

< 

%  =  < mmY  ~  mY)/mY 

and 

(63) 

*»YS  “  VYrJ/vrr,- 

Also,  as  will  be  expressed  by  Eq.  (C6), 

(64) 

x  =  (mi/o-yy0)  -  OOO3), 

and 

(65) 

A 

K  -  0(10).. 

(66) 

When  no  target  is  in  the  sensor’s  field  of  view,  "ideal"  receiver  performance  is  expressed  in 
terms  of  the  estimation  errors  as 


l«*J  «  jt'-ouo-3) 


!%l  «  -  OHO'1). 


From  Zqs.  (62),  (67),  and  (68) 


u(i)  =  K. 
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Thus  Eqs.  (67)  and  (68)  assure  that  u(t)  is  rendered  time  invariant  by  the  AT  processor  when 
no  target  is  in  the  sensor’s  field  of  view. 

From  Eqs.  (41)  ard  (69), 

ntj  -  FAR  •  TOO  -  mJo  exp(-tf2/2)  •  TU\).  (70) 

Equation  (70)  is  identical  to  the  crossing-count  expression  for  uniform  scenes,  Eq,  (49). 

The  kind  of  processor  just  described  is  called  a  constant-false-alarm-rate  (CFAR)  proces¬ 
sor,  since  the  threshold-crossing  rate  is  now  independent  of  time;  that  is,  a  crossing  is  no  more 
likely  to  occur  during  scanning  of  a  region  of  nonuniform  background  than  during  scanning  of  a 
region  of  uniform  background.  With  reference  to  Fig.  4b  for  example,  the  crossing  rate  is  now 
no  greater  at  tp  than  at  any  other  time. 

Unfortunately  the  CFAR  processor  is  generally  a  nonrealizable  ideal.  Depending  on  the 
structural  features  present  in  the  background  radiance  distribution,  large  estimation  errors  (A> 

and  may  be  unavoidable,  giving  rise  to  appreciable  time  dependence  in  Eq.  (62)  for  «(/), 
We  will  show  numerical  examples  that  illustrate  this  point  in  Figs,  lib,  12,  14,  and  15.. 

When  a  target  is  in  the  sensor’s  field  of  view,  "ideal"  receiver  performance  is  obtained 
when  the  time  interval  /,  bracketing  the  target  scan  contains  a  mean-crossing  time  tmi ,  for 
which  (Eqs.  (42)  and  (39)) 

«(/„,)-  0,  (71) 

As  we  discussed  following  Eq.  (42),  a  mean  c;ossing  assures  mat  [1] 

Pn-mj-  1.  (72) 

as  desired, 

Mean-Video  Estimators 

Comparison  of  Eqs.  (67)  and  (68)  shows  that  the  requirement  on  for  ideal  receiver 
perfe.mance  is  two  orders  of  magnitude  more  sevce  than  the  requirement  on  .  In  this 

sense  the  mean-video  estimator  may  be  deemed  a  more  critical  component  than  the  rms 
quantum-noise  estimator. 

In  the  following  sections  we  will  concentrate  our  discussion  of  mean-video  estimators  on 
three  types  of  linear  filters  (Table  1),  evaluating  their  performance  against  a  succession  of 
increasingly  complex  model  backgrounds: 

•  uniform  scenes  of  known  radiance, 

•  uniform  scenes  of  unknown  radiance, 

•  slowly  varying  nonuniform  scenes,  and 

•  scenes  containing  objects  of  arbitrary  contrast  and  size  (such  as  a  target). 
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Table  1  —  Candidate  mean-video  estimators.  All  of 
the  adaptive-threshold  signal  processors  we  consider  in 
this  report  have  the  generic  form  indicated  in  Fig.  5,  in¬ 
corporating  a  filter  for  estimating  the  expected  value  of 
the  video  signal  K(r),  We  consider  the  use  of  three 
types  of  linear  filters  as  candidate  mean- video  estimr.-. 
tors  in  order  to  illustrate  our  approach  to  processor 
design  and  evaluation. 


Signal  Processor 

Mean-Video  Estimate  my(t) 

Fig.  6 

0 

Fig.  10 

YU~Td) 

Fig.  13 

1/2  [Y(t-Td)  +  Y(t+Td )] 

Greatly  improved  mean-video  estimators  can  no  doubt  be  developed  by  incorporating  actual 
measured  background  scenes  into  the  processor-design-and-evaluation  procedure  and  by 
expanding  the  scope  of  candidate  processor  designs  to  include  nonlinear  algorithms  as  well  as 
linear  ones. 

UNIFORM  SCENES  OF  UNKNOWN  BRIGHTNESS 

la  this  section  we  will  illustrate  the  potential  advantage  of  adaptive-threshold  (AT)  pro¬ 
cessing  over  fixed-threshold  (FT)  procc  sing  by  analyzing  IRST  performance  against  the  sim¬ 
plest  possible  type  of  incompletely  specified  scene,  the  uniform  scene  of  unknown  brightness. 

We  will  first  show  that  the  unknown  brightness  level  can  badly  degrade  the  performance 
of  an  FT  system.  We  will  then  analyze  the  AT  processor  structure  of  Fig.  6  and  show  it  can 
maintain  the  same  false-alarm  rate  (FAR)  against  all  uniform  scenes  regardless  of  brightness. 
In  other  words  the  Fig.  6  processor  is  said  to  operate  in  the  constant-false-alarm-rate  (CFAR) 
regime  against  the  class  of  uniform  backgrounds  of  unknown  brightness. 


Fig  6  —  Candidate  processor  structure  that  achieves  constant-false-alarm-rate 
(CFAR)  performance  against  all  spatially  uniform  scenes  regardless  of  bright¬ 
ness  The  transfer  functions  H(J)  and  H0(J)  are  related  by  Eq  (Bl)  Realiza¬ 
tion  of  requires  either  true  dc  response  or  dc  restoration  Comparing 

this  figure  wi'h  Fig  5  shows  that  the  mean  video  has  been  estimated  as  having 
a  zero  value  A  unique  feature  of  this  processor  is  that  the  expected  value  of 
<>!  is  equal  to  ne  rms  value  of  I  (/)  (as  shown  in  Appendix  B) 
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Fixed-Threshold  Processors 

We  analyze  the  false-alarm-rate  performance  of  an  FT  processor  against  a  uniform  scene 
by  use  of  Eqs.  (49),  (50),  and  (20)  with  the  assumption  that  H( 0)  -  0,  resulting  in 

FAR  -  (A/),WJ  exp  -y  (yJ<rY)2  ,  (73) 

It  follows  from  Eqs.  (26)  and  (A5)  that  the  current  variance  [A2)  in  Eq.  (73)  is  directly  pro¬ 
portional  to  the  scene  radiance  mJW  •  cm-2  ■  sr- 1 J 

(74a) 

where 

C  =  (2e  A/  K0RiAitl),  (74b) 

We  define  the  quantities  K0,  Rh  and  AitX  in  Appendix  A. 

For  illustration  we  assume  that  an  acceptably  low  false-alarm  rate  is  achieved  by  choosing 
the  fixed  threshold  y0  in  Eq.  (73)  as  follows: 

y0  “  8  <r  -  8  (CmLr)in,  (75) 

where  mLC  is  the  radiance  of  the  scene  chosen  to  initially  calibrate  the  IRST  md  a  YC  is  the 
corresponding  value  of  <rY  obtained  from  Eq.  (74),  It  follows  from  Eqs.  (73)  through  (75) 
that  the  false-alarm  rate  is  given  by 

FAR  -  FARC  exp  I32(mt  -  mtr)/mt],  (76) 

where  FARf  is  the  calibration  value  of  FAR:; 

FAR,  s  (A/),*,,  exp(-32),  (77) 

It  follows  from  Eq.  (76)  that  a  20%  increase  in  the  scene  radiance  above  the  calibration 
value, 

mL  -  1.2f»i.o  (78) 

causes  a  200-fold  degradation  (increase)  in  the  false-alarm  rate. 

This  simple  example  clearly  demonstrates  that  large-scale  uncertain. ies  in  the  false-alarm-; 
rate  performance  of  fixed-threshold  processors  are  induced  by  relatively  small  uncertainties  in 
scene  brightness. 

Adaptive-Threshold  Processors 

In  this  subsection  we  analyze  the  false-alarm-rate  (FAR)  performance  of  the  adaptive- 
threshold  processors  depicted  in  Figs,  2  and  6  against  uniform  scenes  of  unknown  brightness. 

As  discussed  in  Part  I  111,  Eqs.  (37)  and  (38)  simplify  as  follows  for  uniform  backgrounds 
for  the  Fig.  2  processing  structure: 
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<rfo0-  «rf(l  +  «)  (79) 

and  op^#»  o-j,  (1  +a3),  (80) 

where  a  is  defined  by  Eq.  (54)  as  the  ratio  of  the  noise  bandwidths  of  H0(f)  and  From 

Eqs.  (35),  (36),  (39),  (48),  (79),  and  (80) 


FAR  -  (A/),mj 


fl  +  a3 

1/2 

exp 

-1 

myo  -  my 

2 

[l+a 

2(1  +  a) 

<Ty 

(81) 


If  //(/)  is  a  bandpass  characteristic,  so  that  H( 0)  -  0,  it  follows  from  Eq.  (20)  that 

my-  tf(0)m, -0.  (82) 


We  assume  that  //0(/)  is  a  low-pass  characteristic,  normalized  analogous  to  Eq.  (24): 

max  //<)(/)  -  Ho(0)  -  1.  (83) 

From  Eq.  (83)  and  by  analogy  to  Eq.  (20) 

myo-  KHo(0)mx  -  Kmx.  (84) 

where  K  is  the  threshold  gain,  as  shown  in  Fig.  2. 


From  Eqs.  (82)  and  (84) 

cry 

From  Eqs.  (85),  (26),  and  (A5) 

my0-mr 

<Ty 


K 2  (,tttX/ (T  y) \ 


2ebf  mx 


“  mL. 


(85) 


(86) 


It  follows  from  Eq.  (86)  that  the  argument  of  the  exponential  in  Eq.  (81)  is  linearly  propor¬ 
tional  to  the  scene  radiance  mL,  Thus  the  false-alarm-ratc  performance  of  the  AT  processor 
depicted  in  Fig.  2,  like  the  FT  processor  analyzed  in  tin  previous  subsection,  strongly  depends 
on  the  scene  brightness.  This  is  due  to  the  incomplete  cancellation  of  the  mx  factors  in  the 
numerator  and  denominator  of  Eq.  (86)  :  my0  is  proportional  to  mx,  but  cry  is  proportional  to 
mx.  We  remedy  this  situation  in  Fig.  6,  by  including  a  square-root  device  in  the  threshold  cir¬ 
cuitry. 


Figure  6  shows  that 

Y0(r)  -  K&yU).  (87) 

From  Eqs.  (60)  and  (87) 

my(t)  -  0  (88) 

for  the  Fig.  6  processor.  It  follows  from  Eqs.  (61),  (82),  and  (88)  that 

u(t)  =-  Hj)).  (89) 
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Importantly  it  is  shown  in  Appendix  B  that  a  Y  as  determined  by  the  Fig.  6  processor  is  an 
unbiased  estimate  for  cr  Y:, 

m^-ory,  (90) 

From  Eqs.  (37)  and  (87) 

<*yy0  “  cr +  X2(<r,r  J° 

From  Eqs.  (66),  (90),  (91),  and  (B6) 

K/^kk0)-1-O(10-'‘) 
for  the  Fig.  6  processor.  From  Eqs.  (89)  and  (92) 

u(t)  ~  K. 

From  Eqs.  (41)  and  (93) 

FAR  -  (A f)rm  exp  (-*2/2)  (94) 

for  the  structure  shown  in  Fig.  6,  where  (A/),ms  is  given  by  Eq,  (48)  and  K  is  the  adaptive- 
threshold  gain  (Fig.  6).  Thus  the  false-alarm  rate  of  the  AT  processor  depicted  in  Fig.  6  is 
independent  of  the  brightness  level  of  the  presumably  spatially  uniform  scene. 

Our  discussion  has  so  far  ignored  a  potentially  important  complication:  Eq.  (35)  is  derived 
in  Part  I  (1]  subject  to  the  requirement  that  YU)  and  y0(r)  ate  both  Gaussian  random 
processes.  However,  although  Y(r)  and  4(f)  in  Fig.  6  are  each  Gaussian,  the  square-root  dev¬ 
ice  ensures  that  y0(f)  is  not  Gaussian.  Thus  Eq,  (35)  is  not  strictly  applicable  to  the  structure 
of  Fig.  6, 

However,  the  following  plausibility  argument  suggests  that  a  rigorous  treatment  of  this 
problem  is  probably  unnecessary.  It  follows  from  Fig.  6  and  Eqs.  (C2)  and  (C6)  that 

Z±  «oucr>), 

myg  mB  2  mA 

However,  from  Fig.  6  and  Eq.  (B3) 

mYn  -  Km-  -  K<ty • 

From  the  preceding  two  equations 

vyjvy  =O(10-3).  (95) 

According  to  Eq.  (95)  the  rms  fluctuation  of  y0(f)  is  three  orders  of  magnitude  smaller  than 
the  rms  fluctuation  of  y(f).  Thus  it  appears  unlikely  that  the  statistics  of  (y  -  y0)  are  appre¬ 
ciably  influenced  by  the  statistics  of  y0.  Consequently  it  also  appears  unlikely  that  a  more 
rigorous  treatment  of  this  problem  would  yield  a  result  significantly  different  than  Eq.  (35). 


(91) 

(92) 

(93) 
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Another  important  point  that  we  brushed  over  is  the  requirement  for  a  low-pass  threshold 
filter  (Eq.  (83)).:  Bandpass  characteristics  (ac  coupling)  eliminate  some  practical  problems,  such 
as  1/f  detector  noise.  However,  the  present  analysis  indicates  that  implementation  of  an 
effective  AT  processor  for  operation  against  uniform  scenes  of  unknown  brightness  requires 
knowledge  of  the  dc  component  of  the  photocurrent.  We  will  find  in  succeeding  sections  of 
this  report  that  dc  coupling  (or  synthetic  dc  restoration  [11,12])  is  also  required  by  IRST  AT 
processors  designed  for  operation  against  structured  backgrounds. 


STRUCTURED  SCENES 

We  have  employed  two  kinds  of  background  models  in  the  analysis  thus  far;  uniform 
backgrounds  of  known  brightness,  and  uniform  backgrounds  of  unknown  brightness.  Moreover 
we  have  obtained  a  processor  structure  (Fig.  6)  that  performs  extremely  well  against  these  types 
of  backgrounds.  In  defining  a  hierarchy  of  increasingly  realistic  (and  complex)  background 
models,  a  reasonable  next  step  is  the  consideration  of  slowly  varying  structured  backgrounds.* 

We  will  first  show  that  the  Fig.  6  processor  performs  poorly  against  the  newly  defined 
class  of  slowly  varying  backgrounds.  We  will  establish  this  point  by  showing  that  the  Fig.  6  pro¬ 
cessor  is  susceptible  to  background-induced  mean  crossings  against  nonpointlike  background 
structures.  The  prediction  of  mean  crossings  is  relatively  simple,  requiring  only  the  calculation 
of  the  my  and  nty0  current  statistics.  Thus  we  will  illustrate  a  general  principle:  point-target 
detectors  that  operate  poorly  are  easily  modeled. 

The  deficiencies  of  the  Fig.  6  processor  suggest  a  modification  in  structure  that  will  lead  to 
another  processor  (to  be  depicted  in  Fig.  10).  Our  analysis  of  the  new  processor  will  show  that 
it  can  suppress  mean  crossings  by  slowly  varying  backgrounds.  Nonetheless  we  will  invoke  the 
concept  of  constant-false-alarm-rate  (CFAR)  performance  to  show  that  considerable  room  for 
improvement  still  remains.  We  will  propose  yet  another  processor  structure  (to  be  depicted  in 
Fig.  13)  toward  the  objective  of  achieving  CFAR  performance  against  slowly  varying  back¬ 
grounds.  The  Fig.  13  processor  employs  a  variant  of  "linear  spatial  filtering"  to  establish  the 
threshold  process  Y0  (t). 

Slowly  Varying  Scenes 

Toward  defining  what  we  mean  by  slowly  varying  scenes,  we  rewrite  Eq.  (27)  in  the  form 
my(t)  -  J  ^  h(X)  mx  (t  -  A)  d\.  (96) 


An  unresolved  point  radiator  in  the  field  of  view  (such  as  a  target)  induces  fluctuations  in 
mx{t)  that  are  as  rapid  as  those  of  the  filter  impulse  response  h(t). t  However,  we  assume 


"The  meaning  of  "slowly  varying"  in  the  present  context  will  be  expressed  quantitatively  in  the  next  subsection  Quali¬ 
tatively  it  means  that  all  objects  in  the  field  of  view  must  be  much  larger  than  a  single  pixel  in  spatial  extent  Whether 
or  not  a  given  scene  is  slowly  varying  thus  depends  on  the  instantaneous  field  of  view  of  the  sensor 
tThis  is  assured  by  the  usual  design  equation  2 1/Td  -  I,  where  Td  is  the  dwell  time  of  the  sensor  and  a/ is  the  noise 
bandwidth  of  //(/)  (101 
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here  that  the  time  variations  in  mxit )  are  slow  compared  to  hit),  that  is,  that  mx  i t )  is  not  tar¬ 
getlike.  This  justifies  writing  a  three-term  Taylor  expansion  for  mxit  -  X)." 


mx(t  -  X)  -  mx(t)  -  kmxit)  +  k2mxit)/2.. 

(97) 

Substituting  Eq.  (97)  into  Eq.  (96)  yields 

mYit)  -  mxit)  /j  -  mxit)  12  +  mxit)I}/ 2, 

(98) 

where  by  definition 

/,  =  «,)  a. 

(99) 

h  =  thit)  dt, 

(100) 

and 

/y-  f*  t2  h  it)  dt. 

»/-oo 

(101) 

It  follows  from  Eqs.  (11)  and  (99)  through  (101)  that 

/,  -  Hi 0), 

(102) 

/j  -  H  i 0)  i-j2nV', 

(103) 

and 

/y  -  -//"( 0)  (27T)~2, 

(104) 

where  the  primes  indicate  differentiation  with  respect  to  electrical  frequency  /.  If  Hif)  is  a 
bandpass  characteristic,  then 

//(0)~0.  (105) 


It  can  also  be  shown  that 


Hi 0)  -  0. 


(106) 


as  a  consequence  of  hit)  being  a  real  function  of  time.  From  Eqs.  (98)  th.ough  (106) 

my(f)  =  -(8irJ)-'  //"(0)  mxit),  (107) 


According  to  Eq.  (107)  the  bandpass  filter  Hif)  double-differentiates  a  slowly  varying  input. 
The  corresponding  approximation  to  Eq.  (29)  for  o-  £  is  easily  shown  to  be 

(Ty(t)-2ebfmx(t)..  (108) 


As  shown  by  Eq.  (B3), 


mitU)  -  tr>  (t). 


(109) 


Thus  from  F.qs.  (108)  and  (109)  and  Fig.  6 

myo  it)  -  Km&>  it)  ~  Ki2e\f)U2  m}/2it).  (110) 


*lt  is  simple  to  include  one  or  more  higher  order  terms  in  Eq  (97),,  however,,  the  three-term  expansion  is  sufficient  to 
illustrate  the  point  that  slowly  varying  mx(ii  can  defeat  the  processor  of  Eig  6  by  inducing  mean  crossings 
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In  the  remainder  of  this  section 

•  We  will  assume  that  the  transfer  function  //(/)  has  a  simple  bandpass  characteristic  (to 
be  given  by  Eq.  (111)). 

•  We  will  assign  the  mean  current  mx(t )  a  waveshape  many  dwell  times  in  duration, 
representative  of  an  extended  nontargetlike  object  in  the  sensor’s  field  of  view. 

•  We  will  evaluate  Eqs.  (107)  and  (110)  for  the  mean  filtered  current  mY(t)  and  mean 
threshold  mYo(t)  respectively  and  will  show  that  mean  crossings  occur  (Eq.  42,  Fig.  3a, 
and  the  accompanying  discussion).:  Each  such  mean  crossing  is  correctly  interpreted  as  a 
false  target  report,  that  is,  as  a  "false  alarm." 

•  We  will  conclude  that  the  Fig.  6  processor  is  not  effective  for  suppressing  false  alarms 
from  nonpointlike  structural  features  in  the  background.  The  only  current  statistics 
required  to  make  this  assessment  are  mY(t)  and  mf(j(i);'  there  is  no  need  to  calculate  the 
other  current  statistics  (o-  Y,  a  Y,  <j  Y(),  cr  y)  or  the  threshold-crossing  statistic  thj.  Thus  the 

processor  is  operating  in  a  basically  nonstatistical  regime:  quantum-noise  effects  (as 
reflected  in  the  magnitude  of  <t  Y  for  example)  are  totally  overshadowed  by  background- 
structure  effects.  This  operating  regime,  designated  as  the  "background-structure-limited" 
(BSL)  regime,  has  the  characteristic  that  each  false  alarm  that  occurs  can  be  associated 
with  a  structural  feature  in  the  background. 


Model  Filters  and  Backgrounds 

We  choose  a  form  of  "raised  cosine"  for  the  model-filter  characteristic  H(f)  (Fig.  7): 


IW)i-  sin2(ir/,/2/).  0 

“  1.  /, 

-  cosjM/-/2)/2(/3-/2)),  h 

-  0.  /, 


/  <  / 1, 

(Ilia) 

/</l. 

(111b) 

/  <  /j. 

(111c) 

(11  Id) 

Fig  7  —  Modified  raised-cosine  model  for  the  filter 
characteristic  //(/)  Ou.'  illustrative  numerical  exam¬ 
ples  in  this  report  all  assume  that  |//f/)|  has  the  form 
depicted  here  Analytic  expressions  for  \H(f)\  and  A / 
appear  as  Eqs  (111)  and  (114).  The  analytic  expres¬ 
sion  for  (A f),mi  (Eq.  (48))  is  more  complicated  than 
Eq  (114)  for  A /  However,  the  numerical  value  of 
(A ./)„„,  is  usually  not  much  different  than  the  vaiue  of 
A  / 
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It  can  be  shown  from  Eqs.  (Ill)  that 


H( 0)  -  H'( 0)  -  0 


H"( 0)  -  tt3/2/3,  (113) 

It  follows  from  Eqs.  (25)  and  (111)  that  the  noise  bandwidth  of  the  model  filter  is 

A/  -  (S/2  +  3/j  -  5/,)/8.  (114) 

An  analytical  expression  has  been  obtained  for  the  rms  bandwidth  defined  by  Eq.  (48)  but  is 
not  given  here. 

Our  model  background,  transformed  from  radiance  to  current  by  Eq.  (A6),  gives  rise  to  a 
mean  detected  signal  of  the  following  parametric  form:. 

mx(t)  -  x0(l  +  xt[l  +  (f/r)3]'1)..  (115) 


The  parameters  x0,  r,  and  X)  correspond  respectively  to  a  constant  background  level,  the 
halfwidth  at  half  maximum  of  an  "object"  in  the  scene,  and  the  object/background  peak 
contrast.  Depending  on  whether  r  is  large  or  small  compared  to  the  system  dwell  time  (and  the 
duration  of  hit)),  Eq.  (115)  for  mx(t)  is  representative  of  a  "clutterlike"  or  a  "targetlike"  object 
respectively.  A  representative  plot  of  Eq.  (115)  is  given  in  Fig.  8.  For  the  present  purposes  we 
assume  r  is  much  longer  than  a  dwell  time,  so  that  mx(t)  is  definitely  not  targetlike. 


Fig.  8  —  Model  nonuniform  background  scene, 
transformed  from  radiance  to  average  current  by  Eq. 
(A6)  Our  illustrative  numerical  examples  in  this 
report  all  assume  that  mx(t )  has  the  form  depicted 
here.  An  analytic  expression  for  mx(i)  appears  as  Eq. 
(115).  The  parameters  x0,  r,  and  xt  correspond 
respectively  to  a  constant  background  level,  the 
halfwidth  at  half  maximum  of  an  'object'  in  the  scene, 
and  the  object/ background  peak  contrast.  Negative- 
contrast  objects  are  modeled  by  choosing  X|  <  0. 
Depending  on  whether  r  is  large  or  small  compared  to 
the  system  dwell  lime,  our  model  mx(i)  is  representa¬ 
tive  of  a  'clutterlike"  or  a  'targetlike'  object  respec¬ 
tively. 


From  Eqs.  (107),  (113),  and  (115) 

mYU)  -  my(0)  [1  -  3(//t)31  [1  +  (f/r)3)'3,  (116) 

where 

mr(0)  -  -W8(/i  t)3-  (117) 

From  Eqs.  (108)  and  (115) 

ay(r)  =  (2 elfx0)u2  (l  +  x,(l  +  (t/r)3|-')l/3,  (118) 

with  A/  being  given  by  Eq.  (114). 
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Figure  6  Processor 

We  see  that  for  the  Fig.  6  processor 

W  -  K  crYU  -  rrf)..  (119) 

Taking  the  expected  value  of  Eq.  (119)  and  using  Eq.  (B3),  we  obtain 

mYoU)  -K  <rY(t  -  Td),  (120) 

where  cr  K  is  given  by  Eq.  (118)  for  the  model  filter  and  model  background  of  the  previous  sub¬ 
section. 


Equations  (116)  and  (120)  are  plotted  in  Fig.  9  with  the  illustrative  parameter  values 


*0  -  10-’  A, 

(121a) 

x,  -  1, 

(121b) 

r  -  12.5  Td  s, 

(121c) 

A/  -  5000  Hz, 

(122a) 

f  \  -  2000  Hz, 

(122b) 

(A f)rm  -  4034  Hz, 

(122c) 

K  -  9. 

(123) 

Fig  9  —  iiiusiraiivc  calculations  of  the  mean  video  «; 
and  the  mean  threshold  mt  for  the  processor  dep¬ 
icted  in  Fig.  6.  Model  background  and  filter  parame¬ 
ters  are  as  specified  by  Eqs  (121)  through  (124)  It  is 
assumed  that  the  sensor  is  scanning  an  object  12  5  tar¬ 
get  widths  in  extent,  which  is  an  object  far  loo  large  to 
be  a  target.  This  figure  shows  that  the  Ftg.  6  proces¬ 
sor  experiences  a  mean  crossing  against  the  nontarget- 
like  object,  giving  rise  to  a  false  alarm.  Thus,  the  Fig. 
6  processor  is  ineffective  against  spatially  structured 
scenes. 


From  Eqs.  (118),  (121),  and  (122)  the  nominal  rms  sensor  noise  is 
cry(oo)  -  (2eA/x0)1/2  -  1.27  x  lO'12  A, 

We  calculate  the  peak  contrast  in  our  model  background  from  Eqs.  (115)  and  (121)  as 

( AW  x  ( 0 )  -  Mt(«>)]  =  .Vy*!  =  1  X  10  9  A. 

Thus  the  peak  contrast  in  the  background  is 

aw v (0>  -  mA oo) 

— - 7-7 - =  790 

(7  >(<») 
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times  larger  than  nominal  noise-equivalent  irradiance  (NEI)  of  the  sensor,  for  the  parameter 
values  we  have  chosen  (Eqs.  (121)  and  (122)),  We  assume  the  system  dwell  time  ii  related  to 
the  noise  bandwidth  A/by  the  usual  search-set  rule  of  thumb  (10) 

Td  -  (2A/)-'  -  l(T4s.  (124) 

According  to  Eqs.  (121)  and  (124)  the  full  width  2r  of  the  object  represented  by  Eq. 
(115)  is  25  dwell  times  in  duration,  and  the  signal  from  an  unresolved  target  is  about  2  dwell 
times  in  duration.  Thus  the  size  of  the  object  represented  by  Eq.  (115)  is  12.5  target  widths. 

As  may  be  seen  directly  from  Fig.  9,  the  extended  object  represented  by  Eqs.  (115), 
(121),  and  (124)  does  give  rise  to  a  mean  crossing.  Furthermore,  as  we  discussed  following  Eq. 
(42),  a  mean  crossing  assures  that  a  threshold  crossing  will  occur.  However,  because  the  object 
under  observation  is  about  12.5  target  widths  in  extent,  this  threshold  crossing  must  be  inter¬ 
preted  as  a  false  alarm. 

In  the  next  subsection  we  show  that  simple  modifications  to  the  structure  of  Fig.  6  result 
in  a  processor  design  that  reliably  discriminates  between  point  objects  and  slowly  varying 
extended  objects,  at  the  expense  of  some  performance  against  uniform  backgrounds. 

Improved  Processors 

Inspection  of  Fig.  9  shows  that  the  mean  threshold  (/)  derived  by  the  Fig.  6  processor 
does  not  increase  quickly  enough  or  strongly  enough  to  prevent  a  mean  crossing.  Preveniion  of 
this  mean  crossing  requires  that  the  mean  threshold  advance  in  synchrony  with  the 

mean  current  mt(r).: 

An  AT  receiver  structure  with  a  more  responsive  threshold  mechanism  is  shown  in  Fig 
10.  Instead  of  Eq.  (120),  one  now  obtains 

m y0(t)  «  »/,(/  -  Tj)  +  AV,(r  -  Tj)  (125) 

We  set  the  delay  time  Td  in  Eq.  (125)  and  Fig  10  equal  lo  the  sensor’s  dwell  time.  We 
presume  that  Eqs.  (116)  through  (118)  for  w,  and  <rt  are  still  valid. 


1 _ J 

Htn 

Y0tt) 


Y  (t) 


Fig  10  -  Candidate  processor  with  a  nonzero  mean-video  estimator  Comparison  with  Fig  5  shows 
that  the  mean-video  estimate  is  simply  a  time-delayed  version  of  the  video  signal  mt(r)  -  F(r  -  Td ) 
The  rms  noise  estimate  is  the  same  as  for  the  Fig  6  processor,  enabling  this  structure  to  achieve  C'FAR 
performance  against  all  unilorm  scenes  regardless  oi  radiance  The  time-delay  element  in  the  threshold 
circuit  is  required  to  prevent  self-thresholding  by  a  target  waveform 
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The  performance  of  the  Fig.  10  processor  is  most  easily  assessed  by  means  of  a  numerical 
example..  Once  again  we  employ  Eqs.  (Ill)  and  (1 15)  as  parametric  models  for 
//(/)  and  mx(t)  respectively  with  parameter  values  specified  by  Eqs.  (121)  through  (124).  The 
corresponding  values  of  mY  and  myQ  are  plotted  as  functions  of  time  in  Fig.  11a. 


Fig.  I  la  —  Illustrative  calculations  of  the  mean 
video  my  and  the  mean  threshold  mY  for  the 

Fig.  10  processor.  Model  background  and  Piter 
parameters  are  the  same  as  for  Fig.  9.  The  sen¬ 
sor  is  again  assumed  to  scan  an  object  N  -  12.5 
target  widths  in  extent.  Comparison  with  Fig  9 
shows  that  a  considerable  performance  improve- 
mem  has  been  achieved,  since  the  mean  crossing 
r  present  in  Fig.  9  has  now  been  eliminated 


Comparison  of  Fig.  11a  with  Fig.  9  shows  that  a  considerable  performance  improvement 
has  been  gained:  the  background-induced  mean  crossing  (the  certain  false  alarm)  present  in 
Fig.  9  has  been  eliminated.  Figure  I  la  shows  that  the  mean  threshold  m >0)  appears  to  accu¬ 
rately  track  the  mean  current  m>(f). 

However,  a  performance  assessment  based  solely  on  visual  inspection  of  Fig.  11a  is 
incomplete  and  misleading.  Once  background-induced  mean  crossings  are  indicated  (as  in  Fig. 
9),  the  analysis  need  go  no  further  :  the  performance  of  a  grossly  inadequate  processor  may  be 
assessed  quite  simply,  without  need  to  evaluate  either  u-yycr  ihj.  However,  the  absence  of 
mean  crossings  (as  in  Fig.  11a)  is  not  necessarily  indicative  of  adequate  performance.  We  show 
this  by  calculating  a  300-fold  increase  in  the  expected  number  of  threshold  crossings  over  the 
time  interval  |//r|  <  2  due  to  the  presence  of  the  nontargellike  object. 

From  Eqs.  (41),  (48),  and  (39) 


nij  =  ( A/) exp[-w’(/)/2]  di 

(126) 

for  slowly  varying  scenes,  where 

u(i)  -  l»i,o(r)  ~  ffi|(l))/(rH|J(') 

(127) 

From  Eq.  (37) 

(128) 

The  quantities  mt(r),  </>(/),  and  (i)  are  calculated  from  Eqs  (116),  (118),  and  (125) 
respectively  for  the  Fig.  10  processor.  Evaluation  of  nijU)  by  means  of  Eqs  (126)  through 
(128)  can  now  be  performed,  once  it  is  shown  how  is  calculated 


From  Eq  (60)  and  Fig  10 

F0(r) -«,(/)  + A*  er,(r)t  (129) 
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where 

•hyU)  -  Y(t  -  Tj) 

(130) 

From  Eq.  (130) 

UAU)  -  < TyU  -  Td). 

(131) 

From  Eqs.  (131)  and  (B3) 

-  tj)- 

(132) 

From  Eqs.  (132)  and  (B6) 

O(10-3)., 

(133) 

From  Eqs.  (129),  (133),  and  (131) 

«ry0(')  *  m<rr(t  -  Td). 

(134) 

Finally  from  Eqs.  (128)  and  (134) 

crlyu(l)  -  arj(t)  +  ir2y(l  -  TJ , 

(135) 

where  tr^r)  and  a>(f  -  Td)  are  calculated  from  Eq  (118). 


The  quantity  u(i )  given  by  Eq.  (127)  corresponding  to  Fig.  11a  is  plotted  as  Fig.  11b.  If 
there  is  no  object  in  the  field  of  view,  so  that  x,  — =  0  in  Eq,  (115),  the  expected  number  of 
false  alarms  during  the  time  interval  |//r|  ^  2  is  given  by 

mj  -  (Af),ms  4t  •  cxp(-K}/4)  -  3.2  x  10  s,  (136) 

where  the  appropriate  values  for  r,  (A and  K  are  obtained  from  Eqs.  (121),  (122),  and 
(123)  respectively.  However,  integrating  Eq,  (126)  numerically  with  i/(r)  as  given  by  Fig.  lib, 
we  find  that 

mj  -  1  x  10  \  (137) 


uv'T 


Fig.  lib  —  Illustrative  calculations  of 
uii)  =  —  mr)/<r  tY0  far  the  same  parame¬ 

ter  va!i,js  as  used  for  Fig.  11a.  The  threshold- 
crossmy  rate  increases  about  four  orders  of  mag¬ 
nitude  in  the  neighborhood  of  the  closest- 
approach  lime  ria  The  false- alarm  probability 
during  a  time  interval  |r/r|  $  2,  as  calculated 
from  Eq.  (126),  shows  a  300-fold  increase  due  to 
the  presence  of  the  nontargetlike  object  Thus 
the  time  variation  in  u(i)  reflects  less-than-ideai 
performance. 


Comparison  of  Eqs.  1 137)  and  (136)  shows  that  the  probability  of  a  false  alarm  has 
increased  more  than  300-fold.  This  is  a  significant  deficiency  in  performance,  since  the  full 
width  of  the  object  in  the  field  of  view  is  12.5  limes  larger  than  a  target. 
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We  examine  the  variation  in  processor  performance  with  object  size  by  defining  the 
parameter  A  as 


A'  —  (t/Tj).  (138) 

Thus  A  is  the  ratio  of  object  size  to  target  size.  In  our  preceding  example  we  assumed  a  value 
A  «=  12.5  (Eq.  ( 1 2 1 c) ) .  Plots  of  »(/)  for  object  sizes  A  “  9  and  A  =  6.5  are  shown  in  Fig.  12. 
All  parameter  values  other  than  t  are  as  specified  by  Eqs.  (121)  through  (124). 

Fig.  12  —  Plot  which  together  with  Fig.  11b 
i  u  Vl  shows  the  effect  of  object  sue  on  the  perfor¬ 

mance  of  the  Fig.  10  processor  All  parameter 
values  except  for  object  size  N  are  the  same  as  in 
Fig.  11b.  An  object  with  size  N  »  9  target 
widths  causes  a  six-order-of-magnitude  increase 
in  the  nominal  false-alarrn  probability  during  the 
time  interval  I  z/r  |  ^  2.  As  we  discussed  in  the 
text,  this  is  established  by  evaluating  Eq.  (126) 
for  mj  with  uU)  as  given  by  the  solid  (A  -  9) 
curve  in  this  figure.  According  to  Eq  (139)  a  6% 
false-alarrn  probability  is  calculated  for  the  9- 
target-width  object.  The  zero  intercept  for 
A/-6.5  indicates  a  mean  crossing.  Hence  the  Fig 
10  processor  suffers  a  certain  false  alarm  against 
an  object  6.5  target  widths  in  extent 


The  A  -  9  curve  in  Fig.  12  shows  even  more  dramatically  than  for  A  =  12.5  that  the 
absence  of  mean  crossings  is  not  indicative  of  adequate  performance.  The  threshold-crossing 
rate,  as  calculated  from  Eq.  (126),  briefly  increases  by  seven  orders  of  magritude  as  the  sensor 
is  scanned  across  the  nonlargetlike  object  9  times  larger  than  a  target  'the  probability  of  a 
quantura-noise-induced  threshold  crossing,  originally  small,  grows  to  significant  amplitude  in 
the  neighborhood  of  the  time  at  which  the  difference  heiwef.i  w.  and  assumes  its 
minimum  value  (the  "time  of  closest  approach"). 

Integrating  Eq.  (126)  numerically  with  u(t)  as  given  by  the  A  =  9  curve  in  Fig.  12,  we 
find 

mj  =  0,06,:  A  -  9,  (139) 

corresponding  to  a  6%  probability  of  false  alarm.  If  theic  were  no  object  in  the  field  of  view,  so 
that  *|  «  0  in  Eq.  (115),  the  expected  number  of  false  alarms  during  the  time  interval 
|  f /r  |  <  2  is  obtained  front  Eq.  (136)  as 

ntj  ~  (9/12.5)  x  3.2  x  10  s  -  2.3  x  10-*..  (140) 

The  factor  9/12.5  in  Eq,  (140)  is  due  to  the  linear  dependence  of  mj  on  r  (Eq.  (136)),  and  the 
linear  variation  of  r  with  A  (Eq.  (138)),  Comparison  of  Eqs.  (139)  and  (140)  shows  that  the 
extended  object,  9  target  widths  wide,  has  caused  a  six-order-of-magnitude  increase  in  the  false 
alarm  probability. 

The  dashed  curve  in  Fig.  12  shows  that  the  Fig.  10  processor  suffers  a  mean  crossing  (arid 
hence  a  certain  false  alarm)  for  A  »  6  5,  that  is,  for  objects  6  5  times  luiger  than  a  target 
Clearly  there  is  still  considerable  room  for  improvement  in  processor  performance. 
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Improved  performance  is  obtained  with  the  processor  structure  depicted  in  Fig.  13,  The 
transfer  function  H0(f)  in  Fig.  13  is  related  to  //(/)  as  follows  (Eq.  (B1 ) ) : 

H0(f)  -  (2A /)-'  //(/)  ®  H(f),  (141) 

corresponding  to  an  impulse  response 

M»)  -  (2A /)-'  hHi),:  (142) 

where  //(f)  and  //0(f)  are  the  Fourier  transforms  of  H{f)  and  H0(f)  respectively. 


Fig  13  —  Candidate  adaptive-threshold  signal-processing  structure  I  he 
estimated  value  of  the  mean  video  m(  is  obtained  as  the  average  of  time-, 
delayed  and  t  .ne-advanccd  versions  of  the  video,,  as  we  will  give  a«  Eq 
(143).  Since  the  auanuty  m>  (/)  -  Y(t)  is  proportional  to  the  discrete-lime 
second  derivative  of  K(r),  this  processor  is  a  Laplacian  linear  ..-ouai  niter, 
modified  to  incorporate  the  "efficient  shot  noise  estimator”  developed  in 
Appendix  B 


From  inspection  of  Fig  13  and  Eq.  (142)  the  mean  and  variance  estimates  are  obtained  as 
V')  -  IN*  -  Tj)  +  YU  +  Ttl))/2  (143) 


a;(f)  =  Wr’(f)  ®U(f  -  TJ  +  X(i  +  Tu))/ 2. 
It  can  be  shown  from  Eqs.  (143)  and  (144)  that 


!>0 (f)  -  “(mf  +  +  *  y(°'?+  +  O’)  ) 


+  y(orj  +  try*). 
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where 


(T  y(t  ±  Tj) 

(142a) 

My  (t  ±  Tj). 

(147b) 

Again  we  examine  the  performance  of  our  latest  candidate  processor  structure  by  consid¬ 
ering  some  numerical  examples.  If  the  model  filter  and  model  background  are  again  assumed 
to  be  "iven  by  Eqs  (111)  and  (115)  respectively,  mY(t)  and  <rY(t)  are  still  given  by  "'ns  (116) 
and  (118).  Equations  (126),  (127),  and  (128)  for  m},  «(/),  and  o-yyj f)  are  also  still  applica¬ 
ble.  However,  Eqs.  (145)  and  (146)  now  replace  Eqs.  (125)  and  (135)  respectively. 


The  first  set  of  numerical  calculations  for  the  Fig.  13  processor  is  shown  in  Fig.  14,  where 
all  parameter  values  except  for  r  are  as  specified  by  Eqs.  (121)  through  (124).  Values  of  r  for 
Fig.  14  are  given  by 

T  =  ST,,  .  /V  -  5.  6.5,  9.  (148) 

where  Tlt  =  10“4  s,  as  given  by  Eq.  (124).  When  nij  is  calculated  by  means  of  Eq.  (126),  the 
result  for  the  expected  number  of  false  alarms  during  the  interval  |r/r|  ^  2  is 

ntj  -  9.4  x  10~\  jV  -  9,  (149) 

corresponding  to  the  solid  curve  ( N  *■  9)  in  Fig.  14. 


Fig.  14  —  Performance  assessment  of  the  Fig.  13 
processor  against  objects  of  various  sizes.  All 
parameters  except  for  object  size  are  the  same  as 
for  Figs.  11  and  12.  The  processor  experiences 
no  difficulty  in  suppressing  false  alarms  for  N  - 
9  (against  the  9-target-width  object),  retains  mar¬ 
ginal  capability  for  JV—  6.5,  and  experiences  a  cer¬ 
tain  false  alarm  for  A— 5. 


Comparison  of  Eqs.  (149)  and  (139)  shows  that  the  false-alarm  probability  of  the  Fig.  13 
processor  is  about  five  orders  of  magnitude  smaller  than  the  false-alarm  probability  of  the  Fig. 
10  processor.  This  is  an  important  improvement  in  performance.  Further  comparison,  of  Figs. 
12  and  14,  shows  that  the  Fig.  10  processor  suffers  a  certain  false  alarm  (a  mean  crossing)  for 
N  —  6.5,  but  the  Fig,  13  processor  does  not.  However,  inspection  of  Fig.  14  shows  that  the 
Fig.  13  processor  suffers  a  certain  false  alarm  for  N  -  5  Uor  an  object  5  times  wider  than  a  tar¬ 
get).  Again  the  latest  processor  (Fig.  13)  has  considerable  room  for  improvement  as  a  detector 
of  unresolved  targets. 
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Noise  Bandwidth  and  Dwell  Time 


In  this  subsection  we  will  show  that  system  performance  against  nonuniform  scenes  can 
be  improved  at  the  expense  of  degraded  performance  against  uniform  scenes.  We  presume  that 
the  signal  processor  is  structured  as  in  Fig.  13. 

The  key  to  improving  this  processor’s  performance  against  extended  (nontargetlike) 
objects  is  to  modify  the  conventional  search-set  rule  of  thumb  relating  noise  bandwidth  and 
dwell  time  [10]  as  given  by  Eq.  (124): 

A/-  (2Td)~l.  (150) 

We  now  consider  the  effect  on  system  performance  of  choosing  a  value  for  A/ four  times  larger 
than  specified  by  Eq.  (150): 

A/  —  M(2Td)~\  AT  —  4.  (149) 

We  refer  to  the  processor  described  by  Eq.  (150)  and  Fig.  13  as  an  "AT— 1  processor"  or  "nar¬ 
rowband  processor,"  and  we  refer  to  the  processor  described  by  Eq.  (151)  and  Fig.  13  as  an 
"AT=4  processor"  or  "wide-band  processor."  We  evaluate  the  relative  performance  of  the  two 
processors  by  extending  the  numerical  example  of  the  previous  subsection. 


We  calculate  Fig.  15  by  means  of  the  same  equations  as  Fig.  14,  except  that  we  employ 
the  following  parameter  values  (descriptive  of  the  A/— 4  processor)  instead  of  Eqs.  M21) 
through  (124): 


and 

From  Eqs.  (151)  and  (153) 


X 

o 

1 

© 

1 

> 

(152a) 

X\  -  1, 

(1 52b) 

t  -  NT,,.  N  -  2.5,  5. 

(152c) 

A /  -  20  kHz, 

(153a) 

f\  -  8  kHz, 

(153b) 

(A /)„„,  -  16.136  kHz, 

(153c) 

K  -  9. 

(154) 

T„  -  lO"4  s; 

(155) 

th?'.  is,  the  sensor  dwell  time  is  the  same  as  for  the  previous  example  (Eq.  (124)). 


Comparison  of  Figs.  14  and  15  shows  that  the  AT=  1  processor  suffers  a  certain  false  alarm 
while  the  performance  of  the  AT=4  processor  is  barely  perturbed  for  an  object  5  target  widths  in 
size  (/V=  5).  Further  inspection  of  Fig.  15  shows  that  the  AT- 4  processor  retains  some  capabil¬ 
ity  even  against  extended  objects  just  2.5  target  widths  in  size.*  increasing  the  noise  bandwidth 
still  further,  by  choosing  AT  >  4  in  Eq.  (151),  results  in  even  greater  selectivity  against 
extended  targets. 

*Wc  cannot  evaluate  sensor  performance  using  the  approximate  Eqs.  (116)  and  (118)  for  smaller  values  of  N  due  to 
our  slowly- varying-background  approximation,  Eq.  (98).  This  limitation  is  due  to  our  desire  to  maintain  computational 
simplicity  in  our  numerical  examples  and  is  not  fundamental  to  the  theory. 
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Fig.  IS  —  Performance  assessment  of  the  Fig.  13 
processor  employing  a  wide-bandwidth  elec¬ 
trical  filter.  The  noise  bandwidth  for  this  figure  is 
4  times  larger  than  for  Fig.  14.  Processor  perfor¬ 
mance  is  barely  perturbed  by  the  presence  of  the 
same  5-targel-width  object  ( /V—  5)  thul  would 
cause  the  narrow-bandwidth  (A/“  1)  processor  to 
experience  a  false  alarm  (Fig.  14).  The  wide- 
bandwidth  (A/- 4)  processor  retains  marginal  ca¬ 
pability  against  an  2.5  target-width  object. 
The  improvement  in  performance  against  nonuni¬ 
form  scenes  for  the  A/** 4  processor  compared  to 
the  A/- 1  processor  is  gained  at  the  expense  of 
degraded  performance  against  uniform  scenes. 


We  now  show  that  the  improvement  in  performance  against  nonuniform  scenes  for  values 
M  >  1  is  gained  at  the  expense  of  degraded  performance  against  uniform  scenes. 

With  the  assumption  that  //(/)  in  Fig.  13  is  a  bandpass  characteristic,  it  follows  from  Eqs. 
(82)  and  (147)  that  for  uniform  backgrounds 

my  mf  »  0,  (156a) 

and 

(T  y  —  <Ty  =  (T  y.  (156h) 

Thus,  from  Eqs.  (145)  and  (156) 

niyo  —  K(T  y  (157) 

against  unstructured  scenes  (for  which  x\  —  0  in  Eq.  (115)).  Also,  from  Eqs.  (118)  and  (115) 
(with  x\  ~  0),  we  find 

try  -  (2ex0A/),/2.  (158) 

Thus,  from  Eqs.  (157),  (158),  and  (151) 

my0  «  K(bf)U2  «  K  MU1.  (159) 

According  to  Eq.  (159)  the  threshold  value  myQ  for  the  A/— 4  processor  is  twice  as  large  as  the 

mean  threshold  for  the  A/- 1  processor.  Consequently  the  peak  target  signal  strength  must  be 
twice  as  great  to  induce  a  mean  crossing  (a  certain  target  detection)  in  the  A/— 4  processor  than 
in  the  M—  1  processor.* 

‘The  structure  of  the  Fig.  13  processor  ensures  that  a  target  signal  does  not  experience  significant  self-thresholding. 
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Alternatively  it  follows  from  Eq.  (159)  that  we  can  compensate  for  a  4-fold  increase  in  A/ 
by  halving  the  threshold  gain  K,  Instead  of  Eq.  (154),  the  equations 

K  -  9/2,  M  -  4,  (160) 

-  9,  M  -  4, 

assure  that  the  A/—  4  and  A/-1  processors  have  equal  values  of  mYo  and  hence  equal  values  of 
target  signal  strength  required  to  induce  a  certain  target  detection. 

From  Eqs.  (8),  (126),  (127),  (128),  (145),  (146),  and  (156) 

FAR  oc  exp(-K2/4)  (161) 

for  the  Fig.  13  processor,  for  operation  against  uniform  scenes.  It  follows  from  Eqs.  ( 1 60)  and 
(161)  that  the  reduced  value  of  threshold  gain  required  to  equalize  target-detection  probabilities 
results  in  a  greatly  degraded  uniform-background  fatse-alarm  rate  for  the  M~  4  processor  as 
compared  with  the  M~  1  processor. 

Analogous  to  our  discussion  following  Eq,  (52)  the  performance  disadvantage  described 
can  be  interpreted  as  a  uniform-background  false-alarm-penalty  (FAP)  that  accrues  to  the  use 
of  wfde-bandwidth  processors  (M  >  1).  For  the  Fig.  13  processor  we  can  show  that 

FAP(A/)  -  1.086  A2(l  -  M~])  dB.  (162) 

For  the  previous  numerical  example  (K  «  9,  M  -  4)  we  find 

FAP  -  66  dB. 

Table  2  summarizes  the  results  of  our  numerical  examples. 


The  Promise  of  Focal-Plane-Array  Technology 

The  performance  disadvantage  incurred  by  quadrupling  A  /  (as  w->  pointed  out  following 
Eq.  (161))  can  be  completely  offset  by  replacing  each  detector  element  with  four  series  TD1 
detectors  (as  determined  by  equations  similar  to  Eqs.  (58)  and  (59)  in  the  subsection  “Time 
Delay  and  Integration"!  Thus,,  improvements  m  foea'-plane-array  (FPA>  technology  (!“!]  'lead¬ 
ing  to  large  numbers  of  TD1  elements  may  afTurd  us  "excess"  uniform-background  sensiti vit> 
that  we  can  intentionally  sacrifice  in  the  interest  of  better  performance  against  structured  back¬ 
grounds. 
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Table  2  —  Summary  of  performance  calculations  of  candidate 
processors.  In  performing  example  calculations,  we  assume  that  the 
filter  characteristic  and  background  variation  have  the  parametric  forms 
indicated  in  Eqs.  (Ill)  and  (115)  respectively..  All  calculations  assume 
numerical  values  for  the  parameters  x0>  xb  K,  and  Td  given  by  Eqs. 
(121),  (123),  and  (124),  Calculations  for  narrow-bandwidth  electrical 
filters,  Eqs.  (122),  correspond  to  the  bandwidth  parameter  M r -  !.  and 
calculations  for  wide  bandwidth  filters  correspond  to  Eqs.  053)  and  M 
*  4.  The  parameter  N  is  the  ratio  of  object  size  to  target  size.  Object 
contrast  relative  to  sensor  noise-equivalent  irradiance  (NEI)  is  equal  to 
790  for  the  1  calculations,  and  395  for  the  Af-4  calculations.  The 
minimum  target  signal  strength  required  to  assure  target  detection  is  9 
times  the  sensor  NEI  for  all  calculations.  The  last  entry  in  the  table 
(N  -  2.5,  M  -  4)  indicates  that  a  clutter  object  2.5  target  widths  in 
extent  and  of  amplitude  395/9  -  44  times  greater  than  the  minimum 
target  contrast  required  for  assured  detection  gives  rise  to  a  12%  proba¬ 
bility  of  false  alarm. 


Processor 


Object  Size  Bandwidth  False*A!arm  Probability,  \th\  <  2 


Nominal 


Degraded 


Fig.  10 


Fig.  13 


5.2  x  10~17  1 

3.2  x  10 -8  lx  10~5 

2.3  x  10~8  6  x  10"2 

1.7  x  10~8  1 

2.3  x  10"8  9.4  x  10“ 7 

1.7  x  10~8  3.5  x  10' 2 

1.3  x  10“8  1 


Fig.  13 


5 

2.5 


5.2  x  10- 8 
2.6  x  IQ’8 


Appendix  A 

CALCULATING  THE  AVERAGE  CURRENT  FROM  BACKGROUND  DATA 


Our  objective  in  this  appendix  is  to  derive  Eq.  (A5),  which  expresses  the  average  value  of 
the  random  current  X(t)  (Fig.  1) 

mx(t)  s  £{*(*))  (Al) 

as  a  function  of  the  scene  radiance  mL  and  a  number  of  important  senior  parameters.  We  will 
derive  Eq.  (A5)  subject  to  the  assumptions  that  the  scene  is  spatially  uniform  (that  the  radiance 
mL  is  independent  of  view  angle)  and  that  the  average  current  mx  is  constant  with  time.  We 
will  then  generalize  Eq.  (A5)  to  Eq.  (A6)  to  allow  for  the  possibility  of  spatially  nonuniform 
background  radiances  and  time- variable  average  currents  mx(t ).. 

The  average  value  of  the  photocurrent  X(t)  (F  g.  1)  may  be  calculated  as 

mx  ~  pe  mp,  (A2) 

where  mx  is  defined  by  Eq.  (Al),  tj  is  the  quantum  efficiency  of  the  detector  (electrons  per 
photon],  e  is  the  electronic  charge  [coulombs  per  electron],  and  is  the  average  photon  flux 
incident  on  the  detector  (photons  per  second].  The  quantities  mx  and  mQ  in  Eq.  (A2)  are  both 
averages  over  the  photon  fluctuation  statistics  of  the  incident  light.  The  meaning  of  the  ensem¬ 
ble  averaging  process  in  the  present  context  is  discussed  in  Part  I  [1];  in  particular,  the  average 
in  Eq.  (Al)  is  not  a  time  average. 

The  mean  photon  flux  incident  on  the  detector,  in  Eq.  (A2),  is  given  by 

m0  “  ^dei  tnE/ln>,  (A3) 

where  is  the  active  area  of  the  photodetector  (cm2],  mE  is  the  focal-plane  irradiance 
iW-cm-2),  and  hv  in  the  radiant  energy  [J]  per  photon. 

The  focal-plane  irradiance  ttiE  in  Eq.  (A3)  is  calculated  as 

mE  =  7rT0mi/(2/#)2  =  K0mL,  (A4) 

where  t0  is  the  transmittance  of  the  optics  (dimensionless],  /#  is  the  focal-length  ratio  of  the 
optics  Idimensionless],  and  mL  is  the  radiance  of  the  scene  (W  cm"2  sr_1l. 

It  follows  from  Eqs.  (A2)  through  (A4)  that 

=  K0  R,  Aiel  m,  (A5) 

where  the  current  responsivity  Rt  is  defined  a-,  =  r)e/hv  iA-W"1].  Equation  (AS)  is  valid  for 
only  uniform  backgrounds  that  fill  the  sensor’s  field  of  view..  Nonuniform  backgrounds  can  be 
characterized  in  terms  of  a  radiance  distribution 

m i  (0)  =  nii  («..  )3).. 
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where  a  and  /3  are  azimuth  and  elevation  angles  respectively.  The  appropriate  generalization  of 
Eq,  (A5)  when  the  background  is  nonuniform  is  [1] 

mx(t)  =  K0R,  f  MTF{i)?(f)m/f)e('2r'  yl)  (if,,  (A6) 


where 


f  =  two-dimensional  spatial  frequency  {(cycles  •  cnTl)2], 

MTF(f)  =  modulation  transfer  function  of  the  optical  train  [dimensionless], 
v  =  focal-plane  scan  velocity  [cm  s  ']  =  2nd/Tr , 


d  =  focal  length  of  the  optics  [cm], 
TF  =  system  frame  time  [s], 


and  F(f)  and  m,(f)  are  the  two-dimensional  Fourier  transforms  of  the  quantities  P( r)  and 
m/f-t/d): 

f(f)  =  J*  P(t)  cxpijlnf  •  r)  dr  (A7) 

and 

my(f)  =  J*  mi  (-t/d)  expijlnf  ■  r)  dr,,  (A8) 

where  r  is  the  two-dimensional  focal-plane  location  vector  [cm].  For  the  detector  geometry 
depicted  in  Fig.  Al,  the  quantity  P(r)  in  Eq.  (A7)  is  defined  as 

P(t)  «  1,  rE  (A9a) 

-  0 ,  r  4  (A9b) 

It  follows  from  Eqs  (A7)  and  (A9)  that 

iF(0)  =  /tdc,  [cm2].  (A10) 

For  uniform  backgrounds 


Fig.  Al  -  Focal-plane  geomelry  The 
focal-plane  irradiance  is  sianonary  in  the 
coordinate  system  with  origin  0  Vector 
r„  locales  the  center  of  a  detector  of  area 
For  scanning  sensors  r„  is  a  func-. 
tion  ol  lime 
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my( f)  -  mt  8(f),,  (All) 

where  8(-)  is  the  Dirac  delta  function.  From  Eqs  (A6),  (A10),  (All),  and  the  conventional 
MTF  normalization, 


it  follows  that 


MTF(O)  =  lt 


(A12) 


mx  *  A'„  R,  /f 0cl  mt  (A13) 

for  uniform  backgrounds.  Equation  (A13)  is  consistent  with  Eq,  (A5),  as  expected.  Finally, 
we  note  that  the  quantities  m,  (0)  and  f),  like  m*(r),  are  defined  as  ensemble  averages  over 
the  photon  fluctuation  statistics  of  the  incident  light. 


3d 


Appendix  B 

EFFICIENT  ESTIMATORS  FOR  NONSTATIONARY  SHOT  NOISE 
INTRODUCTION 

A  main  result  of  this  report  is  the  proposal  of  a  new  signal  processor  structure  for  generat¬ 
ing  high-confidence  estimates  of  the  time-varying  rms  quantum  noise  induced  in  a 
background-limited  (BLIP)  sensor  as  it  scans  a  nonuniform  scene.  This  new  signai-processing 
structure  has  the  general  form  illustrated  in  Fig.  Bl.  For  reasons  that  we  will  explain  shortly, 
the  processor  depicted  in  Fig.  Bl  is  called  an  "efficient  shot-noise1'  (ESN)  estimator. 


Fig.  Bl  —  Efficient  estimator  for  nonstationary  shot  noise.  The 
transfer  functions  H(f)  and  H0(/)  are  related  by  Eq.  (Bl) 
Realization  of  Ho (/)  requires  either  true  dc  response  or  dc  res¬ 
toration  (Eq,  B2).  An  important  feature  of  this  processor  is  that 
the  expected  value  of  a  Y  is  equal  to  the  rms  value  of  T(r)  (Eq. 
B3).  All  of  the  candidate  processors  in  this  report  use  the  struc¬ 
ture  shown  here  to  perform  the  rms  quantum-noise  estimation 
function  indicated  in  Fig  5. 


The  quantity  (>,(r)  in  Fig  Bl  is  an  estimate  for  the  rms  value  of  the  random  process 
>'(();  that  is,  rr  > (/)  is  an  estimate  for  trY(i).  Inspection  of  Fig.  Bl  shows  that  i >t  is  esta¬ 
blished  by  cascading  a  linear  filter  //„(/)  with  a  square-root  device  and  an  amplifier  of  gain 
(2 elf)112.  The  frequency  characteristic  //„(/)  is  given  by 

//„(/)  -  (2A/)  1  //(/)  «//(/),;  (Bl) 

where  ®is  the  convolution  operator  (Eq.  (10)).  Also,  e  is  the  electronic  charge  (1.6  x  10  19  C) 
and  A/ is  the  noise  bandwidth  of  #(/),  as  defined  by  Eq.  (25).  The  transfer  function  H0(f) 
given  by  Eq.  (Bl)  is  always  a  low-pass  function,  From  Eqs.  (25)  and  (Bl) 

Ho(0)  -  max  H0(f)  =  1,  (B2) 

/ 

Thus,  realization  of  Wy(/)  always  requires  a  sensor  with  either  true  dc  response  or  dc  restora¬ 
tion  [11,12]. 
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The  noise  estimate  fry  established  by  the  Fig.  B1  processor,  like  X '(f)  and  Y(t),  is  a  ran¬ 
dom  process.  We  show  in  this  appendix  that 

It!-  (t)  =  f  { CT  >  ( /  )  }  =  <Ty(t),:  (B3) 

"  y 


where  £{•)  denotes  an  ensemble  average  over  the  photon-fluctuation  statistics  of  the  incident 
light  [1],  Equation  (B3)  indicates  that  the  ensemble  average  value  of  the  noise  estimate  fry  is 
equal  to  the  actual  rms  noise  o-  that  is,  that  fr  y  is  an  "unbiased"estimate  for  o-  Y .  Even  so,  fr  t 
might  assume  a  sufficiently  large  spread  of  values  about  its  average  to  make  fry  a  poor  estimate 
for  cry.  To  ensure  that  d>  is  not  a  poor  estimate,  we  calculate  the  mean-square  estimation 
error: 

(t\ ^  =  E[(fr  y  -  (T  y)1]-  (B4) 

We  show  in  this  appendix  that 

“  f2  A/  Vo  <B5) 

for  the  Fig.  B1  processor,  where  e  is  the  electronic  charge,  A/is  the  noise  bandwidth  of  //(/), 
and  A/o  is  the  noise  bandwidth  of  H0(f)  (Eq.  (BD). 

It  follows  from  Eq.  (B3)  that  fr  y(t)  is  generally  a  nonstationary  random  process,  since  the 
mean  value  of  fr  Y  is  generally  nonconstant  with  time.  It  is  thus  surprising  to  see  fiom  Eq.  (B5) 
that  the  variance  of  fry  is  time  independent. 

It  follows  from  Eqs.  (B5)  and  (B3)  and  from  looking  ahead  to  Eqs.  (C4)  through  (£6) 
that  the  relative  mean-square  estimation  error  of  the  Fig.  B1  processor  is  extremely  small: 

“  0(10-*),  (B6) 

According  to  Eq.  (B6)  the  variance  of  fr  Y  is  six  orders  of  magnitude  smaller  than  the  variance 

of  F(f).  This  results  in  significant  simplifications  in  some  of  the  analysis  (as  we  will  show  fol¬ 
lowing  Eq.  (95)  for  example.) 

One  commonly  used  approach  to  evaluating  the  effectiveness  of  an  estimator  invokes  the 
Cramer-Rao  inequality,  which  provides  a  lower  bound  on  o-j(  [4],  We  evaluate  the  Cramer-Rao 

lower  bound  for  in  this  appendix  and  show  it  to  equal  crj(  as  given  by  Eq.  (B5).  Con¬ 
sistent  with  the  conventional  terminology  of  statistics,  it  follows  that  the  Fig.  B1  processor 
establishes  an  "efficient"  estimate  for  a  y. 

Before  proceding  with  the  statistical  arguments  described,  we  will  first  show  that  the 
discrete-time  implementation  of  Fig.  B1  takes  the  particularly  simple  form  shown  in  Fig  B2. 
The  quantity  7>  in  Fig.  B2  is  both  the  dweil  time  of  the  sensor  and  the  integration/sample  time 
of  the  postdetector  filter. 

DISCRETE-TIME  IMPLEMENTATION  OF  THE  ESN  ESTIMATOR 

We  now  assume  that  the  postdetector  filter  is  simply  an  integrate,  sample,  and  hold  dev¬ 
ice  The  waveform  at  the  output  of  the  postdetector  filter  is  the  discrete-time  random  process 
MU),  where 
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Fig.  B2  -  Discrete-time  implementation  of  the  efficient  shot-noise  estima¬ 
tor.  The  generic  structure  depicted  in  Fig.  B1  simplifies  considerably  when 
it  is  implemented  as  a  discrete-time  device,  since  we  now  have 
//(/)-  The  quantity  Td  is  both  the  dwell  time  of  the  sensor  and 

the  integration/ sample  time  of  the  postdetector  filter. 


W(k) 


T- 1 
•  d 


X 


(*+l/2)T. 


X(t)  dr. 


(B7) 


With  the  function  red  (■)  defined  as  shown  in  Fig.  B3  [13),  Eq.  (B7)  can  be  written  as 

W(k)  =  HkTj), 

where 

Hi)  -  Tj'  red  (t/Tj)  «>X0). 

, RECT(X) 


(B8) 

(B9) 


Fig,  B3  -  Definition  of  the  function  reel  (*) 

The  impulse  response  of  the  postdetector  filter  is 
a  scaled  version  of  reel  (*)  when  the  Fig.  B1  pro¬ 
cessor  is  implemented  in  discrete  time  (Eq  (B9)). 

Comparing  Eqs.  (10)  and  (B9),  we  see  that  the  impulse  response  of  the  postdetector  filter  is 

hit)  =  Tj1  red  UlTj).  (BIO) 

Taking  the  Fourier  transform  of  Eq.  (BIO),  we  obtain  [13) 

H(f)  =  sine  (JTd)  =  (ir/T,,)-1  sin  (nfTd).  (BID 
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It  follows  from  Eq.  (Bll)  that  H(f)  is  normalized  properly  (Eq.  (24)): 

tf(G)  -max  //(/)-  1,  (B12) 

f 

Thus,  from  Eqs.  (21),  (25),  and  (BIO), 

A/-(2rrf)-',  (B13) 

In  addition  it  follows  from  Eqs.  (142),  (BIO),  and  (B13)  that 

h0U)  -  hit).  (B14) 

Finally  it  follows  from  Eqs.  (B13)  and  (B14)  that  the  discrete-time  implementation  of  Fig.  B1 
takes  the  particularly  simple  form  depicted  in  Fig.  B2. 

ESTIMATOR  BIAS 

Our  purpose  in  this  section  is  to  show  that 

£|<r  y  —  <r  rl  “  0.  (B15) 

According  to  Eq.  (B15)  the  mean  estimation  error  is  zero,  so  that  fry  provides  an  unbiased  esti¬ 
mate  for  a  r •• 

It  follows  from  Fig.  B1  that 

frY(t)  -  (2eif)'n  Bit),  (B16) 

Taking  the  expected  value  of  both  sides  of  Eq.  (B16),  we  obtain 

mi>(  -  (2eA/)l/2  mB.  (B17) 

From  Eq.  (B17)  and  the  relation  -  m\n  (1  -«),  which  we  will  give  as  Eq.  (C40),  where  e 
-  (1/2)  ( <rA/2mA )2,  which  we  will  give  as  Eq.  (C42), 

-  {2e A/  mA)U2  (1  -  «).•  (B18) 

It  follows  from  inspection  of  Fig.  B1  that 

A(t)  -  Mr)  ©.K(r).:  (B19) 

where  ®is  the  convolution  operator  and  //0(r)  is  the  Fourier  inverse  of  //u </). 

h0(t)  -  F  'l//o(/)|..  (B20) 

It  can  be  shown  from  Eqs.  (Bl)  and  (B20)  and  the  convolution  theorem  for  Fourier  transforms 
that 

/i„</)  -  J7A/1  1  i:Ht).  (B21) 

Thus  from  Eqs  (B19)  and  (B21) 

A(i)  =  (2 A/)  1  h2(i)  ®X(t).  (B22) 

Taking  the  expected  value  of  both  sides  of  Eq.  (B22),  we  find 

/«,(/)  =  (2 \J)  1  h2(t)  (B23) 
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Substituting  Eq.  (B 23)  into  Eq.  (B18),  we  obtain 


(t)  **  I eli2(i)  ®mx(t)]'n  (1  -  e), 

(B24) 

However,  from  Eq.  (29) 

( Ty(i )  -  eh2(i) 

(B25) 

Thus  from  Eqs.  (B24)  and  (B25) 

m ,-r>,(t)  "  (Ty(t) (1  -  «),, 

(B26) 

or 

vr-n.r, 

- -  -  t  -  O(10~6). 

O'! 

(B27) 

According  to  Eq.  (B27)  the  mean  estimation  error  for  cry  is  six  orders  of  magnitude  smaller 
than  cry  itself.  The  Fig.  B1  processor  is  thus  described  as  an  asymptotically  unbiased  estimator 
for  a  y. 

MEAN-SQUARE  ESTIMATION  ERROR 


In  deriving  Eq.  (B5)  we  first  note  that 

=  £(((>,  -  <r  t)J|  -  £{<r>)  +  cr\  -  2<r  >  £((>,) 

It  follows  from  Eq.  (B26)  that 

£{<>rl  =  m&)  -  cr > (1  -  J, 
where  again.,  as  we  will  give  as  Eq.  (C42), 

«  =  y  Ury/2m4)\ 

Since  from  Fig.  B1 

a  ](i)  =  2c M 

it  follows  from  Eq.  (B22)  that 

<?;(/)  =  elr(t)®X(t). 

From  Eq.  (B32) 

£{«>;(/)}  =■  <•  lr(i)  ®  £{.V(t)} 

=  e  lr(t)  ®m \(t). 

From  Eqs  (B33)  and  (B25) 

£{(> ;  (i )  |  =  i r  j  (r ). 

From  Eqs  (B28),  (B29),  and  (B34) 

<r ;  =  2  e  <r j . 


(B28) 

(B29) 

(B30> 

(B31) 

(B32) 

(B33) 

( B  3  4 ) 

(B35) 
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Since,  as  we  will  give  as  Eq.  (C6),  aAlmA  -O(10~3),  it  follows  from  Eqs.  (B35)  and  (B30) 
that 

(°"ff  ~  O(10-6),,  (B36) 

verifying  Eq.  (B6). 

From  Eqs.  (B30),  (B23),  (B25),  and  equations  we  will  derive  as  Eqs.  (C40)  through 
(C43) 

2«  -  e2  A/  A/o/cr  j.,  (B37) 

Equation  (B5),  which  is  one  of  the  principal  results  of  this  appendix,  now  follows  from  Eqs. 
(B35)  and  (B37):] 


*lt  -  e2  A/  A/0. 


(B38) 


CRAMER-RAO  INEQUALITY 


As  shown  by  Van  Trees  [4,  p.  66),  the  mean-square  estimation  error  defined  by  Eq.  (B4) 
must  satisfy  the  Cramdr-Rao  inequality 


.2  > 

.i  ^  a  r. 


CRi.H- 


(B39) 


The  Cramdr-Rao  lower  bound  <rlRlB  is  defined  as 

=  -E{d2r.  ,,  >n/dr|,r(ayky)l.  (B40) 

where  dvt.„r  denotes  the  second  partial  derivative  with  respect  to  ay  and  /(.  |  (cr  >  | cr is 

the  probability  density  function  for  conditioned  on  rr,.  Since  <r(  is  obtained  as  the  square 
root  of  a  Gaussian  random  process,  its  density  function  is  readily  adapted  from  the  iiterature 
(14,  pp.  287-290] 


o {  -  try 


2^,  o, 


+  <t> 


-  2  . 

a)  + 

Oj 

2(r 

"t 

o> 

(B41 ) 


where  is  given  by  Eq.  (B38),/*(  )  is  the  unit-step  function,  and  0(  )  is  defined  as 

4>(x)  "  (2rr)~l/2  exp(~x2/2).  (B42> 

Since 


+  a] 

2<r 

-,<r> 

* 


2  »»■„ 


(B43) 


for  all  value  of  a)s  it  follows  from  Eq  (B6)  \  lat  the  second  term  within  the  brackets  on  the 
right-hand  "ue  of  Eq.  (B41)  contributes  negligibly  >o  fCj  ,  so  that 


a  2  -  a  y 


2  a 


(B44) 
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Taking  the  natural  logarithm  of  Eq,  (B44)  and  then  the  second  partial  derivative  with  respect  to 
cry,  we  obtain 

-rty..r  On  fiy k>>  -  tl  +  n&y/cry)'}.  (B45) 

Also  it  can  be  shown  from  Eq.  (B44)  that 

E[{aYfav)A\-\,  (B46) 

to  the  same  order  of  approximation  as  in  Eq.  (C27), 

From  Eqs.  (B40),  (B45),  and  (B46) 

*Irlb  ”  vly-  A/  A/0.  (B47) 

Since  the  mean-square  estimation  error  o-?  is  equal  to  the  Cramer-Rao  lower  bound,  the  Fig. 
B1  processor  is  a  generator  of  "efficient"  estimates  for  a  y. 
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SQUARE  ROOT  OF  A  GAUSSIAN  RANDOM  PROCESS 
INTRODUCTION 

Our  purpose  in  this  appendix  is  to  derive  simple  expressions  fcr  the  mean  and  variance  of 
the  square  root  Bit)  of  a  Gaussian  random  process  A  (i)  (Fig.  6).  The  mean  and  variance  of 
A  (t)  are  called  mAit)  and  j -A(t)  respectively.  The  mean  and  variance  of  Bit)  are  called  mBit) 
andoj(r)  respectively. 

We  will  derive  the  relationships 

mBit)  »  mj/2(r)  (Cl) 

and 

<rB(t)/mB(<)  -  y  [<rA(t)/mA(t)\,  (C2) 

subject  to  the  assumption  that 

tTA(t)/mA(t)  «  1.  (C3) 

The  validity  of  Eq.  (C3)  is  assured  in  that  the  process  Ait)  is  the  output  of  a  background- 
limited  (BLIP)  infrared  sensor,  for  which 

o- \  =  2fA/0  mA(  (C4) 


where  e  is  the  electronic  charge  and  A/a  is  the  noise  bandwidth  (Eq  (25))  of  the  filter  HoU) 
(Fig.  6).  Although  Eq.  (C4)  is  only  approximately  correct  for  time- varying  mA(t)  (as  we  will 
see  from  Fq.  (108)),  it  can  be  used  to  provide  an  order-of-magmiude  eslimate  for  <r  JmA.  We 
assume  for  illustrative  purposes  that 

mA  =  10“ 9  A  (C5a) 


and 


A/o  =  2500  Hz, 


It  follows  from  Eqs.  (C4)  and  (C'5)  that 


justifying  f"q  (C'3). 


(r  4/ m  4 


odo-’),. 


(C5b) 

(C6) 


MATHEM\fICAL  PRELIMINARIES 


We  write  the  probability  density  function  of  Ait)  as 


J4(a)  =  <j 


a  -  m . 


(T  . 


iCV 
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where  we  define  <(>  (■)  as  the  normalized  Gaussian  density: 


< i>(x )  =  (2ir)~i/2  exp  ( ~x2/2).~ 

(C8) 

Since 

*(!)-  U0)|1/2, 

(C9) 

it  follows  from  Eq.  (C7)  that 

\a\U2<f>  da. 

(CIO) 

Also, 

al  -  £{£2}  -  mj, 

(Cl  1) 

where  £{■}  is  the  statistical  expectation  operator  and 

£(fi2}  »=  cr^1  J  \a\4>  da. 

J-oc 

(C12) 

We  devote  the  remainder  of  this  appendix  to  developing  estimates  for  Eqs.  (CIO)  and  (Cl 2). 

EVALUATION  OF  £{£2! 


(C13). 

(C14) 

(Cl  5) 

(C16) 

(C17) 

(C18) 

(Cl  9) 
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We  now  write  Eq.  (Cl 5)  as 

r 00 

h  “  cr^1  J0  a  ^  b/vA^sia)  da, 

where  by  definition 


g(a)  =  </> 


a  +  mA 


<*A 


/ 


It  can  be  shown  that 


max  g(a)  -  g(0)  -  (2tt)i/2  <f>(mA/(rA). 

0<a<» 


From  Eqs.;  (C20)  and  (C22) 

h  <  ^a1  (2 tt)1/2  4>{mAl<rA)  Jo  a<l>(a/(TA)  da. 
Since  the  integral  in  Eq,  (C23)  evaluates  exactly  as  o-j,  it  follows  that 

h  <  (2 tt)1/2  (t a  (f>(mA/<rA). 

From  Eqs.  (Cl 9)  and  (C24) 

£{B2}  -  mA(\  +  tA), 

where 


€  A  <  2(2tt)  1/2 


&A 


mA 


<t> 


f»A 


<TA 


(C20) 

(C21) 

(C22) 

(C2.‘) 

(C24) 

(C25) 

(C26) 


It  follows  from  Eqs.  (C6),  (C2S),  and  (C26)  that  the  quantity  eA  is  truly  negligible  and  that  to 
all  intents 


E{B2)-  mA< 


(C27) 


EVALUATION  OF  m, 


We  now  write  Eq.  (CIO)  as 
where 


mB  -  /3  +  /4, 


and 


h  =  *4l  /o 


a  ~  mA 
<ta 


da, 


a  +  mA 

0‘A 


da. 


(C28) 

(C29) 


(C30) 


Equation  (C30)  can  be  estimated  by  means  similar  to  those  used  in  estimating  Eq.  (CIS)  to 
obtain  the  bound 


/  4< 


mA 


1/2 


'Jtt 

<*A 

2 

mA 

2 

mA 

<*A 

(C31) 
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It  follows  from  Eqs.  (C6)  and  (C31)  that 


h  =  0, 


(C32) 


to  the  same  order  of  approximation  as  Eq.  (C27), 

The  first  step  in  evaluating  13  is  to  write  Eq.  (C29)  as 


;3 


1  + 


o'  A 


mA 


1/2 


(j){u)du. 


Writing  a  binomial  expansion  for  the  square  root  of  the  bracketed  quantity  in  Eq. 
obtain 


h  “  mj 


1/2 


,  1 

a  A 

,  1 

O’A 

2 

ho  +  ~2 

"t A 

/31_  8 

™a 

hi  + 

where 


4  =  >“1-2-3' 

It  can  be  shown  from  Eqs.  (C6)  and  (C35)  that 

ho”"  ^32  =  1 
and 

/3,°o 

to  the  same  order  of  approximation  as  Eqs.  (C27)  and  (C32). 

From  Eqs.  (C28),  (C32),  (C34),  (C36),  and  (C37) 


mB  = 


.  1 

"A 

2 

8 

mA 

From  Eqs.  (Cl  1),  (C27),  and  (C38) 


O’  B  =  mA 


2m  A 


2  o-j 


4mA  ' 


From  Eqs.  (C6),  (C38) ,  and  (C39) 


=  mP  (1  -  «) 


mu  =  mA 


and 


where 


(Tti/mB=  —  ( (rA/mA )  (1  +«).: 


«  =  y (tr 4/2 mA)2=  0(  10  b). 


(C33) 
(C33),  we 

(C34) 

(C35) 

(C36) 

(C37) 


(C38) 


(C39) 


(CO> 


(C41 ) 


(C42) 


Thus,  Eqs.  (Cl)  and  (C2)  have  been  validated,  and  the  relative  errors  have  been  estimated 
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Finally  from  Eqs.  (C4)  and  (C39) 

2  _  mA 

4  mA 


eAfo/2. 


(C43) 


I(  is  interesting  to  compare  Eqs.  (C4)  and  (C43).  According  to  Eq.  (C4)  the  noise  at  the  input 
to  the  square-root  device  is  a  function  of  mA  (and  hence  the  mean  photon  flux  incident  on  the 
detector).  However,  Eq.  (C43)  indicates  that  the  noise  at  the  output  of  the  square-root  device 
is  independent  of  the  mean  photon  flux. 
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